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ABSTRACT 

Though the small-scale, low-ionization knots, filaments and jets (LISs) of planetary 
nebulae (PNe) are known for ^30 yr, some of their observational properties are not 
well established. In consequence our ability to include them in the wider context of the 
formation and evolution of PNe is directly affected. Why most structures have lower 
densities than the PN shells hosting them? Is their intense emission in low-ionization 
lines the key to their main excitation mechanism? Therefore, if considered altogether, 
can LISs line ratios, chemical abundances and kinematics enlighten the interplay be¬ 
tween the different excitation and formation processes? Here we present a spectroscopic 
analysis of hve PNe that possess LISs conhrming that ah nebular components have 
comparable electron temperatures, whereas the electron density is systematically lower 
in LISs than in the surrounding nebula. Chemical abundances of LISs versus other 
PN components do not show signihcant differences as well. By using diagnostic dia¬ 
grams from shock models, we demonstrate that LISs’ main excitation is due to shocks, 
whereas the other components are mainly photo-ionized. We also propose new diag¬ 
nostic diagrams involving a few emission lines ([N ll], [O ill], [S ll]) and \og{ishocks/^ 
where ishocks and C are the ionization photon fluxes due to the shocks and the central 
star ionizing continuum, respectively. A robust relation differentiating the structures 
is found, with the shock-excited clearly having \og{ishocks/^i.)>-^] while the photo- 
ionized show \og{i shock A transition zone, with -2< log{fshocks/h)<-^ where 
both mechanisms are equally important, is also dehned. 

Key words: ISM: abundances - ISM: jets and outflows - ISM: kinematics and 
dynamics - planetary nebulae: individual: NGC 6572, IC 4846, K 1-2, Wray 1-17, 
NGC 6891 


1 INTRODUCTION 

Over the last 30 years, imaging surveys of planetary nebu¬ 
lae (PNe) have brought to light shapes and structures with 
high complexity (Manchado et al. 1996; Gorny et al. 1999; 
Boumis et al. 2003, 2006; Parker et al. 2006; Sahai 2011; 
Sabin et al. 2014). Several works have been carried out thus 
far attempting to group PNe based on their morphologies 
(Balick 1987; Schwarz et al. 1993; Manchado et al. 1996; Sa¬ 
hai et al. 2011). Any morphological classihcation of PNe 
relies on the shape of the large-scale structures (nebular 
components) such as rims, attached shells and halos, and 
divide them into hve main classes: point-symmetric, ellipti¬ 
cal (that includes round), bipolar, multi-polar and irregular. 
These large-scale structures are better identihed in hydro¬ 
gen recombination lines as well as in bright forbidden 
emission lines. 


The formation of spherically symmetric PNe can ad¬ 
equately be explained under the interacting stellar winds 
model proposed by Kwok et al. (1978). Although, the de¬ 
viation from spherical symmetry is still an open question 
(see Balick & Frank 2002; Shaw 2012). A considerable ef¬ 
fort has been made to explain the formation of complex 
PNe morphologies. Most models consider the presence of 
an equatorial density enhancement in order to collimate the 
fast stellar wind and result in the formation of the axis- 
symmetric PNe (e.g. Mellema et al. 1991; Mellema 1995). 
The former equatorial density enhancement is usually at¬ 
tributed to the mass-exchange interaction between the com¬ 
ponents in a close binary system (Soker & Livio 1994; Nord- 
haus & Blackman 2006). Single rotating stars, with or with¬ 
out magnetic helds, have also been investigated as potential 
origin of the equatorial density enhancement and therefore, 
the formation of aspherical PNe (Garcfa-Segura et al. 1999). 
Although, very recent magneto-hydrodynamic models have 
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shown that complex PNe, such as bipolar, cannot be descen¬ 
dant of single stars (Garcia-Segura et ah 2014). 

Another enigmatic and poorly understood component 
in PNe are the structures on smaller scales than the main 
nebular components (Balick et ah 1998; Corradi et ah 1996; 
Gongalves et ah 2001). These structures are seen mostly in 
the light of low-ionization lines, such as [N ll], [S ll], [O ii] and 
[O l]. Gorradi et al. (1996) unveiled their presence in various 
PNe, independently the morphological type, by using the 
image [N il]/[0 ill] line ratio technique, which emphasizes 
the excitation degree through the nebulae. 

These small-scale, low-ionization structures (hereafter 
LISs; Gongalves et al. 2001) exhibit a variety of morpholog¬ 
ical types such as knots, jets, jet-like, filaments and head- 
tails (Balick et al. 1998; Gorradi et al. 1996; Gongalves et 
al. 2001), whilst they cover a wide range of expansion veloc¬ 
ities from few tens to hundreds of km s“^. Based on their 
expansion velocities, they are labelled as fast, low-ionization 
emission regions (FLIERs; Balick et al. 1993), bipolar, rota¬ 
tion episodic jets (BRETs; Lopez et al. 1995) or slow mov¬ 
ing low ionization emitting regions (SLOWERs; Perinotto 
2000). Gongalves et al. (2001), studying a large number of 
PNe with LISs, reached the conclusion that LISs appear in 
all morphological classes of PNe, suggesting that they may 
not be associated with the mechanisms responsible for the 
deviation of PNe structures from the spherical symmetry. 

Several models have been proposed to explain the for¬ 
mation of LISs, but none of them can provide a general 
mechanism for all the different types. Jet formation in PNe 
is usually associated with the interaction of a binary sys¬ 
tem and an accretion disk around a companion (e.g. Soker 
& Livio 1994; Blackman et al. 2001). As for the knots, any 
clue for their formation also needs to consider if they appear 
in pairs or are randomly distributed (Gongalves et al. 2001). 
Symmetric pairs of low-ionization knots seem to appear due 
to dynamical or/and radiation instabilities (Garcfa-Segura 
et al. 1999). Eor the isolated knots, in situ instabilities Soker 
& Reveg (1998), and the stagnation knot’s model (Stef¬ 
fen et al. 2001) seem to be plausible formation mechanisms 
that explain some of their characteristics, as much as the 
shock models from Dopita (1997). Eor more details on for¬ 
mation mechanisms of LISs see the comprehensive review 
by Gongalves et al. (2004). A common characteristic of the 
models for the formation of LISs is the determination of their 
age with respect to the main body of the nebula. It has been 
observationally found that some jets predate the main nebu¬ 
lar structure, whereas others appear to be younger, like that 
of NGG 6337 (Tocknell et al. 2014; Jones et al. 2014) and the 
jets and knots of the IG 4846, Wray 17-1 and K 1-2 studied 
in this paper. Altogether the latter two papers quoted, and 
previous results for the age of LIS in PNe (Gongalves 2004) 
reveal a likely connection of LISs with the evolution of the 
PNe central stars. 

Miszalski et al. (2009a) had also proposed a possible 
link between LISs and the central stars of PNe. These au¬ 
thors argue that LISs seem to be very common in PNe with 
hydrogen-deficient stars, e.g [WR]. In this scenario, LISs are 
likely associated with the turbulence of the strong stellar 
winds from these, though this connection is still lacking fur¬ 
ther confirmation. In addition, the same authors claim that 
there may also exist a relation between the presence of LISs 


and the binarity of the central stars, with 40% of post com¬ 
mon envelope (GE) PN possessing LISs. 

The physical and chemical properties e.g. electron tem¬ 
perature (Te), electron density (Ae), ionic and elemental 
abundances of LISs have also been studied by different 
groups over the past e decades (Balick et al. 1993, 1998; 
Hajian et al. 1997; Gongalves et al. 2003, 2004, 2009; Leal- 
Eerreira et. al. 2011; Monteiro et al. 2013). No significant 
difference between the Te of LISs and that of the higher 
excitation nebular component (e.g. core, rims, shells) has 
been found, contrary to the Ae, which is found to be sys¬ 
tematically lower in LISs when compared to the surrounding 
medium. This result is in conflict with the theoretical pre¬ 
dictions from the formation models of LISs, which predict 
LISs to have several times higher density than the surround¬ 
ing ionized medium. Gongalves et al. (2009) claim that the 
mass of LISs must be mostly neutral in order to overcome 
the discrepancy between models and observations, since all 
the formation models of LISs calculate the total density of 
gas (dust, atomic and molecular) and not the electron den¬ 
sity, which corresponds only to the ionized fraction of the 
gas. 

The strong [N ll] emission line found in LISs compared 
to the surrounding ionized nebula was firstly attributed to 
a significant local overabundance of nitrogen (Balick et al. 
1994). Mellema et al. (1998) and, more recently, Gongalves 
et al. (2006) reached the conclusion that an overabundance 
of N is not a necessary condition for getting such a strong 
line-emission. Based on photo-ionization models, Gongalves 
et al. (2006) found that the N+/N = 0^/0 assumption of 
the ionization correction factor (IGE, Kingsburgh & Barlow 
1994) should be wrong, and resulting in seemingly overabun¬ 
dance of N in LISs. 

Accordingly to the theoretical work by Aleman & Gru- 
enwald (2011), the peak intensity of the low-ionization op¬ 
tical lines ([O ii], [N ii], [S ii], [N i] and [O i]), as well as the 
H 2 co-rot at ional lines, occurs in the same narrow transition 
zone between the ionized and neutral (photo-dissociation) 
regions. Given that molecular hydrogen can be excited ei¬ 
ther by shocks or by absorbing UV-photons emitted from 
the central star, the enhancement of low-ionization lines 
in LISs may be attributed to these mechanisms. Raga et 
al. (2008) show that shock models with low and moderate 
photo-ionization rates are able to reproduce emission line ra¬ 
tios typical of shock-excitation regions, whereas those mod¬ 
els with high photo-ionization rates result in emission line 
ratios similar to photo-ionized nebulae. Therefore, besides 
the expansion velocities, the evolutionary stage (e.g. the age) 
and the photo-ionization rate (or the strength) of the central 
star are crucial parameters in order to distinguish between 
predominantly shock- and photo-excited structures. 

In this series of papers, new spectroscopic data from a 
sample of 10 PNe that possess LISs are combined with kine¬ 
matic data in order to study, at the same time, the physico¬ 
chemical and morpho-kinematic properties of LISs and other 
nebular components. The first 5 PNe, with LISs embedded 
in the main nebular shells, are presented here. The remain¬ 
ing PNe will appear in a forthcoming publication. The ob¬ 
servations and data analysis are described in Section 2. The 
emission line fluxes, the physico-chemical properties, and the 
morpho-kinematic characteristics of the PNe themselves, as 
well as of their LISs, are presented in Section 3 and 4, re- 
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wavelength (A) 

Figure 1. Observed spectra of the NW jet in Wray 17-1 (upper 
panel), covering the wavelength range of 3600-8000 A, obtained 
with DFOSC@Danish, and SE outer NEB in NGC 6891 (lower 
panel) with a wavelength range of 3650-7050 A, and observed 
with the IDS attached to the Isaac Newton telescope. 


spectively. The results of this work are discussed in Section 
5, and we wrap up with the conclusions in Section 6. 


2 OBSERVATIONS 

2.1 Intermediate-dispersion Spectroscopy 

Intermediate-dispersion spectra of three PNe (NGC 6891 
and NGC 6572, IC 4846) were obtained with the 2.5 m Isaac 
Newton Telescope (INT) at the Observatorio del Roque de 
los Muchachos, Spain, on 2001 August 30, 31 and September 
4, respectively. The Intermediate Dispersion Spectrograph 
(IDS) was used with the 235 mm camera and the R300V 
grating resulting in a scale of 3.3 A pixel”^ and covering the 
wavelengths range of 3650-7000 A. The spatial scale of the 
instrument was 0.70 arcsec pixel”^, with the TEK5 CCD. 
The slit width was 1.5 arcsec and the slit length 4 arcmin. 

The spectra of K 1-2 and Wray 17-1 were obtained 
with the 1.54 m Danish telescope at the European South¬ 
ern Observatory (ESC), La Silla (Chile), on 1997 April 10 
and 11, respectively. The Danish Faint Object Spectrograph 
and Camera (DFOSC) was used with the 2000 x 2000 CD 
(15/im pixel) and the Grism #4 (300 lines mm”^), resulting 
in 2.2 A pixel”from 3600 to 8000 A. The spatial scale of 
the instrument was 0.40 arcsec pixel”^. The slit width and 
length were 1.0 arcsec and >13.7 arcmin, respectively. 

Spectro-photometric standard stars were also observed 
with the INT and Danish telescopes, in order to calibrate the 
spectra, while the fluxes were corrected for the atmospheric 
extinction and interstellar reddening. During the night, bias 
frames, twilight and tungsten flat-held exposures, and wave- 
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length calibrations were also obtained. Data reduction was 
performed using the standard IRAF instructions for long-slit 
spectra. Individual images were bias subtracted and flat-held 
corrected using a series of twilight flat frames. Two exam¬ 
ples of the spectra under analysis, one for each telescope 
configuration, are presented below in Fig. 1. 


2.2 High-dispersion Spectroscopy 

High-dispersion, long-slit spectroscopy in the Ho+[N ll] 
and/or [O ill] emission lines, for 4 out of the 5 PNe in our 
sample, were acquired from the SPM Kinematic Catalogue 
of Galactic Planetary Nebulae (Lopez et al. 2012; hereafter 
KCGPN). The spectra of three PNe (NGC 6572, NGC 6891 
and K 1-2) were obtained at the 2.1-m, f/7.5 telescope in 
San Pedro Martir National Observatory (Mexico), whereas 
for Wray 17-1, the data were obtained from the 3.9-m, f/8 
Anglo-Australian telescope in the Siding Spring Observa¬ 
tory (Australia), and in both cases by using the Manchester 
Echelle Spectrometer (MES-SPM; Meaburn et al. 2003). 
Two slit widths of 70 and 150/im were used, correspond¬ 
ing to a velocity resolution of 9.2 and 11.5 km s”^, respec¬ 
tively. Additional kinematic information were also drawn for 
the literature. 


3 RESULTS: PHYSICO-CHEMICAL 
PROPERTIES 

In the even numbered tables, from Table 2 to 16, we list 
the emission line fluxes, corrected for atmospheric extinction 
and interstellar reddening, for several nebular components, 
e.g. LISs (knots, jets, filaments), inner and outer main nebu¬ 
lar regions (NEB) and for the entire PN, whenever possible. 
All fluxes are normalized to F(H/3)=100. The interstellar 
extinction c(H/3) was derived from the Balmer Ha/H/3 ratio 
(eq. 1), using the interstellar extinction law by Fitzpatrick 
(1999) and Rv=3.1, 




1 lo >(H«)/-F(H/?) 

0.348 ^ 2.85 


( 1 ) 


where 0.348 is the relative logarithmic extinction coefficient 
for H/3/Ha. 

Te and for each nebular component were derived by 
using the temden task in iraf (Shaw & Dufour 1995) and 
are presented in the lowermost part of the flux tables. Given 
that the abundances derived from collisionally excited lines 
are strongly dependent on Te, two combinations of Te and 
Ne have been used depending on the degree of ionization 
of each ion, whenever possible. The combination of Te[N ll] 
and Ae[S ii] covers the low-excitation regions (N^, and 
S^), whereas the combination of Te[0 III] and Ae[S ll] cov¬ 
ers the medium-excitation regions of these PNe (e.g. 

Ne^^). Significant differences between the two diagnos¬ 
tic temperatures are found in IC 4846, K 1-2 and NGC 6572, 
which make the usage of the two diagnostic temperatures 
necessary for properly deriving the ionic and chemical abun¬ 
dances. 

In addition, for a proper estimation of Ne from the 
[Ar iv] diagnostic lines, we calculated and subtracted the 
contribution of the He I A4712 recombination line and those 
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of the [Ne iv] AA4724, 4726 emission lines to the [Ar iv] 
A4711 one. In particular, the contribution of the He i A4712 
line was calculated using the theoretical work by Benjamin 
et ah (1999), whereas the contribution of the [Ne iv] AA4724, 
4726 emission lines were considered negligible in all cases, 
given that the central stars of the PNe in this work are not 
hot enough to ionize Ne^^ (P.L=97.11 eV), except the cen¬ 
tral star of Wray 17-1 (Teff=140000 K; Rauch and Werner 
1997). Although, the of this nebula is extremely low, 
close to the lower limit, and small changes in the [Ar iv] 
4711/4741 line ratio due to the contribution of the [Ne iv] 
lines will not have a significant impact in the Nq. In short, 
a comparison of derived from the different diagnostics 
lines does not show significant differences that could be as¬ 
sociated with the [Ne iv] lines. We have to mention here the 
extreme case of NGC 6891, for which we find significantly 
higher Ae[Ar iv], compared to Ae[S ll] and Ae[Cl III], by 
factors of 6 and 3.5, respectively. This probably indicates 
strong density stratification, with denser inner regions with 
respect to the outer ones. 

3.1 IC 4846 

IC 4846 (PN G027.6-09.6) is a relatively compact nebula 
with a size of 10 arcsec. Low angular resolution images of 
this nebula (Miranda et al. 2001) show two ellipsoidal shells 
at PAs of 11° (inner shell) and 54° (outer shell). These au¬ 
thors also show the [N Iij/Ha line ratio image which reveals 
the presence of features enhanced in N^: i) a pair of knots, 
attached to the inner shell, at PA=11°; and ii) a pair of 
filaments, attached to the outer shell, at PA=54°. 

A higher resolution Ha image of IG 4846 (Fig. 2), ob¬ 
tained with the HST in 2001 (program ID:8345; PI: Sahai), 
clarifies that the pair of filaments is made of highly colli¬ 
mated structures (jets) with some indications of precession. 
In this new image, the pair of knots along PA=11° appears 
to be the result of the interaction between the inner and 
outer shell. It is worth mentioning here the astonishing sim¬ 
ilarity between IG 4846 and Fleming 1 (Boffin et al. 2012). 
Both PNe exhibit a pair of knotty, processing, bipolar jets 
or BRETS (Lopez et al. 1995), likely associated with an ac¬ 
cretion disk around a binary system. This binary system has 
already been confirmed in Fleming 1 (Boffin et al. 2012) but 
not yet in IG 4846. 

Nevertheless, the central star of IG 4846 has been clas¬ 
sified as a weak emission line or WEL-type star (Acker et 
al. 1992; Gorny et al. 2009). The typical line features of 
WELs and/or [WG] stars, such as the A4650 (G III-G IV) 
and the A5805 (G IV) lines, were detected in our spectra 
with a full width at half maximum (EWHM) of 22±2 and 
24±2 A, respectively. These detections, therefore, reinforce 
the possibility the central star of IG 4846 may belong to 
the rare group of binary systems with a WR or WEL-type 
companion. 

In Table 2 we list the emission line fluxes of the inner 
shell, the pair of jets at PA=54°, and the integrated emission 
of the nebula along the slit (entire nebula). Unfortunately, 
no spectra along PA=15° were obtained. The exact loca¬ 
tion of the extracted portions of the spectra are shown, and 
labelled, in Eig. 2. 

IG 4846’s interstellar extinction, c(H/3), is found to vary 
between 0.48±0.03 and 0.53±0.05, in excellent agreement 


Ha 



offset (arcsec) 


Figure 2. Upper panel: The HST Ha image of IC 4846. The size 
of the field shown is 10x10 arcsec^. The nebular components un¬ 
der analysis are indicated in the image. The extraction windows 
for these components are indicated by the corresponding boxes: 
the entire nebula (entire NEB: 10 arcsec, in red), the inner neb¬ 
ula (inner NEB: 3.5 arcsec, in cyan) and the NE and SW jets 
(2.1 arcsec, in blue). North is up, and East to the left. Lower 
panels: The radial profile of the [N ll] and [O ill] emission-line 
fluxes, normalized to 1.0 and the [N ii]/[0 ill] line ratio, along 
the slit. 


with those derived from previous studies (0.47, Kaler 1986; 
0.50, Gahn et al. 1992). Notice, however, that the most re¬ 
cently published value, by Wesson et al. (2005), is higher by 
a factor of 1.45. One possible explanation for this discrep¬ 
ancy may be the different slit positions (PA=0° in Wesson 
et al. 2005). The logarithmic H/3 fluxes are calculated equal 
to -11.62, -13.49 and -13.45 and 11.54 for the inner shell, NE 
jet, SW jet, and the entire nebula, respectively. Acker et al. 
(1991) give two H/3 fluxes for this nebula (-11.56 and -11.26), 
which were derived from the Haute-Provence Observatory 
spectroscopic data (OHP-GGD) with a slit width of 2.5 arc¬ 
sec, and from the ESO-GGD, spectra with a slit width of 
4 arcsec, respectively. Our value for the entire nebula bet- 
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Figure 3. Upper panel: Ho; image of Wray 17-1 obtained from 
the ESO Archive. The field shown is 110x110 arcsec^. The neb¬ 
ular components under analysis are indicated in the image. The 
extraction windows for these components are indicated by the cor¬ 
responding boxes: the inner nebular regions (inner NEB: 11.6 arc¬ 
sec, in red); the NW and SE blobs (11.6 arcsec and 10.8 arcsec, 
in blue); the knots (NW knot of 4.4 arcsec, SE knot of 3.6 arcsec, 
NE knot of 2.8 arcsec, and SW knot of 3.6 arcsec, in green); and 
the arcs (E arc of 14.8 arcsec and W arc of 6 arcsec, in cyan). 
North is up, and East is to the left. Lower panels: The radial 
profile of the [N ll] and [O ill] emission-line fluxes, normalized to 
1.0 and the [N ll]/ [O ill] line ratio along the slits. Panels (a) and 
(b) correspond to PA=155°, whereas (c) and (d) to PA=75°). 


ter agrees with the value derived from the OHP-CCD data 
because of the same slit widths. 

The lower panels in Fig. 2 show the fluxes, normalized to 
1.0, of the [N il] (solid-line) and [O ill] (dashed-line) emission 
lines along the slit (panel a), as well as the [N ll]/[0 III] 
line ratio (panel b). The (spatial) line profiles and the ratio 
of these two lines are important in order to better define 
the sizes of each structure along the slit. The size of each 
nebular component is indicated in the figures by the coloured 
regions. 

Ne and Te of IC 4846 are estimated using the diagnostic 
line ratios of sulfur, argon, chlorine and oxygen. No signifi¬ 
cant difference is found among them (see Table 2). The com¬ 
parison of Ne and Te between the nebular components shows 
that the jets are significantly less dense than the inner shells, 
by a factor of 4, whereas, within the errors, all structures 
have the same Te. Our Ne values are slightly higher than 
those found by Wang et al. (2004). We have to notice here 
that an extremely high Ae[S il] of ~20000 cm“^ has been 
reported by Wesson et al. (2005), which is not confirmed by 
our analysis, neither for the jets nor for the inner shell. Our 
analysis reveals that the jets in IC 4846 have low Ae, with 
respect to the main nebular component (inner NEB), and 
therefore effaces any argument that the LISs in this nebula 
could have Ne higher than the surrounding medium. The 
questionable line ratios given by Wesson et al. (2005) were 
also previously mentioned by Delgado-Inglada et al. (2015). 

In Table 3 we present the ionic abundances (obtained 
using the IONIC task in IRAF; Shaw & Dufour 1995) as well as 
the total chemical abundances of IC 4846 for each nebular 
component. The total ones were computed using the ion¬ 
ization correction factor (ICE) from Kingsburgh & Barlow 
(1994). These authors do not provide an ICE for Cl, thus we 
used the equation by given Liu et al. (2000). 

The low He abundance and N/0 ratio of IC 4846 indi¬ 
cate a non-type I nebula. We do not find any trend in chem¬ 
ical abundances among the nebular component; therefore, 
they are all constant within the uncertainties. A compari¬ 
son between our total chemical abundances and those found 
in the literature show a good agreement within the errors, 
except for Ar (Perinotto et al. 2004; Wesson et al. 2005). 
The low N abundance in the SW jet is highly uncertain due 
to the ICE scheme for N (mentioned in the Introduction), 
as well as to the fact that only N^ was measured in the 
different regions. However, both N/H of the SW and NE jet 
agree, within the errors. 


3.2 Wray 17-1 

Wray 17-1, or PN G258.0-15.7, displays a diffuse complex 
nebula of approximately 80 arcsec in diameter (Fig. 3). A 
pair of bright regions is seen in the [O ill], [N ll] and [S ll] 
emission line images embedded in the diffuse nebula, along 
PA=155°. The [N ll]/[0 ill] line ratio image of the nebula 
unveils highly collimated jet-like structures (Corradi et al. 
1999). At the end of these jet-like structures (hereafter fil¬ 
aments), a pair of knots can also be discerned, whereas a 
second pair of knots is also apparent along the PA=75°. 
The latter is found to be surrounded by two much fainter 
arc-like structures. These arcs show an enhanced [O III] /Ha 
ratio, suggesting a shock interaction (Guerrero et al. 2013). 
The central star of Wray 17-1 has been found to be a H- 
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deficient P Cygni star with a Teff=140000 K and log(g)=6.3 
(Rauch & Werner 1997). The formation of jets and/or knots 
may be associated with the strong turbulent stellar winds 
from this star (Miszalski et ah 2009b) or a putative accre¬ 
tion disk around a close binary system (Soker & Livio 1994, 
Miszalski et ah 2009b). 

Several nebular regions/components are selected for the 
analysis of this nebula: the pair of filaments, the two pairs 
of knots, the two arcs and four inner nebular regions, for 
direct comparison with the properties of the entire nebula. 
In Fig. 3, we label each of these features with the spatial re¬ 
gions under analysis. The lower panels show the line profiles 
of the normalized fluxes of the [N il] and [O ill] emission 
lines, as well as their ratio. The position of the filaments 
and knots is noticeable in these diagrams. The emission line 
fluxes, the absolute flux of H/3, c(Hyd), Ne and Te for all 
these regions are presented in Tables 4 (PA=155°) and 6 
(PA=75°). The c(Hyd) of this nebula varies from 0.04 to 
0.13 among the nebular component. The average value is 
0.08±0.03, which is in a good agreement with the value de¬ 
rived by Tylenda et al. (1992). Te [0 III] and Te[N il] are 
found to vary from 11100 K to 14700 K and from 11625 K 
to 13625 K, respectively. Within the uncertainties, all these 
measurements agree with one another. 

Unlike the Te, where no variation among the nebular 
components is found, may show some variation. In par¬ 
ticular, the LISs (knots and jet) have higher Nq compared 
to the inner nebular regions, but because of their high un¬ 
certainties, we can say whether the LISs are, in general, less 
or more dense than the inner nebular regions. Moreover, the 
Ne values of the inner nebular regions are close to the low 
density limit of 100 cm“^, and they should be taken with 
caution. 

We present here the first measurements of the ionic and 
total chemical abundances of Wray 17-1 (Tables 5 and 7). 
Surprisingly, the knots of Wray 17-1 do not present any 
He^ (all He is doubly ionized) and hence the abundance 
of He cannot be calculated properly, similarly to K 1-2 (see 
below). The high effective temperature of the central star 
(Te//=140000 K) is consistent with the high excitation of 
the nebula. Unlike the knots, the filaments in Wray 17-1 
appear to have much lower excitation, where He^ is more 
abundant than He^^. 

3.3 K 1-2 

K 1-2 is an elliptical nebula with very low surface brightness 
and a known binary system. Highly collimated structures 
with several knots on opposite directions of the central star 
and a jet-like structure oriented at PA=153° can easily be 
discerned in the [N ll] and [S ll], as well as the [O III] and 
hydrogen recombination line images (Corradi et al. 1999). 
These authors also found several other knots almost per¬ 
pendicular to the direction of the jet-like structure, similar 
to the case of Wray 17-1. The post-CE close binary system 
found at the centre of K 1-2 would be more easily consistent 
with the presence of the collimated jet-like structures and 
knots, as argued in the theoretical works from Soker and 
Livio (1994), a number of other later works, and the recent 
observational works by Tocknell et al. (2014) and Jones et 
al. (2014). 

Our spectroscopic data of K 1-2 were obtained only 
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Figure 4. Upper panel: [N ll] image of K 1-2 obtained from 
the ESO archive. The size of the field shown is 62x53 arcsec^. 
The nebular components under analysis are indicated in the im¬ 
age. The extraction windows for these components are indicated 
by the corresponding boxes: the Northwest nebula (NW NEB: 
35.6 arcsec, in red), the NW filament (7.6 arcsec, in blue) and the 
knots (NW and SE knots: 3.6 arcsec, in green). North is up, and 
East is to the left. Lower panels: The radial profile of the [N ll] 
and [O III] emission-line fluxes, normalized to 1.0, and the [N ll]/ 
[O III] line ratio, along the slit. 


along the PA=153° with the results presented in Table 8. 
The nebular regions under analysis in this nebula are 
shown in Fig. 4. We calculate c(H/3) between 0.22±0.06 and 
0.26±0.06 for the various nebular components (NW jet, NW 
knot and SE knot), which are close to the value derived by 
Exter el al. (2003) for both north and south jets, namely 
0.25. It is worth mentioning here that the c(H/3) of the en¬ 
tire nebula has been estimated to be 0.11 (Exter et al. 2003), 
which suggests that the jet and knots are dusty features, at 
least much more than the surrounding nebula. 

The absolute H/3 fluxes for each nebular component are 
also given in Table 8. For the NW NEB component, the 
T(Hyd) is 11.4 xl0“^^ ergs s“^ cm“^, which is in very good 
agreement with the value of the northern jet derived by Ex¬ 
ter et al. (2003). We also confirm that Te[0 ill] is ~5000 K 
higher than Te[N ll], as previously found, but our values are 
2000 K and 1000 K, respectively, lower than those in Exter 
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et al. (2003). Ne is found to vary between 422zbll0 cm“^ 
and 888zbl95 cm“^ and, due to the large uncertainties, no 
unambiguous evidence of variation in Ne between the LISs 
and nebulae can be provided. 

The lower panels of Fig. 4 displays the radial prohles 
of the [N li] and [O ill] emission-line normalized fluxes, as 
well as the [N ii]/[0 iil] line ratio, along the slit. The exact 
position of the knots can easily be discerned. Moreover, one 
can also see that the NW jet is a separate component from 
the NW knot. The [N il]/[0 ill] is found to be very high at 
the position of the knot, close to 0.8, and decreases outwards. 
The positions of two knots (small green regions) with offsets 
of 8 and 14 arcsec, respectively, are barely seen in the [N ll] 
line, but they are further away with a sharp peak in the 
[N ii]/[0 III] line ratio. These two knots have been noticed 
by Corradi et al. (1999). 

Our NW NEB is the same region studied by Exter et 
al. (2003), and the chemical abundances are found to be in 
a good agreement, in the two works. NW NEB has high 
He++/H+ and low 0+/0++ ionic abundance ratios indicat¬ 
ing a high excitation nebula. Surprisingly, both knots, NW 
and SE, do not show He^ lines. Moreover, the high excita¬ 
tion He^^ and ions are found to be more abundant in 

these components. Why these knots are so highly excited is 
an open question. A possible explanation is that the contri¬ 
bution of the main nebula in these components is significant, 
especially in high-ionization lines (Exter et al. 2003). 

3.4 NGC 6891 

NGC 6891 (PN G054.1-12.1) displays a bright inner prolate 
elliptical shell with the major axis oriented along PA=135°, 
and an inclination angle of 50-55°, surrounded by a diffuse 
spherically symmetric nebula (Guerrero et al. 2000; Palen et 
al. 2002). Guerrero et al. (2000) pointed out the presence of 
a (collimated) outflow between the inner and outer nebulae, 
as well as a pair of faint knots at the tip of these outflows, 
due to the interaction with the outer nebula (see Eig. 1 from 
Guerrero et al. 2000). However, HST images of this nebula 
do not reveal the presence of any collimated outflow (Eig. 5). 
Erom our analysis, we found these regions as normal outer 
nebular regions (outer NEB), instead of LISs. 

The central star of NGG 6891 has been classihed as an 
0-type star (McGarthy et al. 1990; Guerrero & De Marco 
2013) and as a WEL star (weak emission line, Tylenda et 
al. 1993). The strong P Gygni prohle found in its spectrum 
implies a terminal velocity of 1200-1400 km s“^(Marcolino 
et al. 2007). Neither the P Gygni nor any wind variability 
have been detected in the central star spectrum, by Guer¬ 
rero & De Marco (2013). Nevertheless, the detection in its 
spectrum of the A4650 (G III 4647/51 and G IV 4658) and 
A5805 (G IV 5801-5812) line features, in conjunction with 
the absence of the emission line G III A5696 line, support the 
WEL-type classihcation. The He ii AA4541, 5412 absorption 
lines, previously reported by Parthasarathy et al. (1998) and 
Marcolino V de Araujo (2003), were also detected by us. 

Two spectra were obtained for this nebula, with the 
slits positioned along PA 45° and 135°. The identihcation 
of all nebular components used for the analysis is given in 
Eig. 5. The logarithmic E(H/3) flux derived for the entire 
NEB are -11.29 and -11.31 for the PA=135° and PA=45°, 
respectively, in good agreement with those from Acker et 


al. (1991) work. c(H/3) is calculated between 0.03±0.02 and 
0.06±0.04 (for the different nebular components). Be aware 
that, our values are substantially lower than those derived by 
Kaler et al. (1970) and Wesson et al. (2005), namely 0.41 and 
0.29, respectively. The different slit positions between the 
aforementioned studies and ours may cause this discrepancy 
in c(H/3). 

Te[0 III] of the entire nebula, for both PAs, are 
9600±1630 K and 9500±1420 K, whereas Te[N ii] is found 
to be higher by -3000 K (12720±6360 and 12190±5900 K). 
These value agree with the previously published ones: 
Te[0 III] =9200 K and Te[N il]=10600 K (Henry et al. 
2004) and Te[0 iii]=9330 K (Wesson et al. 2005). As for 
the Ae, we derive Ae[S ii] = 1290±390 cm“^, Ae[Gl iii] = 
2460±490 cm“^, and Ae[Ar iv]=2508±640 cm“^, which also 
agree within the errors with those densities derived by Wes¬ 
son et al. (2005). The high Ae value of 10000 cm“^, reported 
by Henry et al. (2004), is not confirmed in this work, nor in 
any other so far available in the literature (Kaler 1970; Wes¬ 
son et al. 2005). 

A comparison of Te and Ae among the inner NEB, outer 
NEB, and knots, do not reveal any important difference. The 
surprisingly high Ae[Ar iv] found in the inner NEB, may be 
indicating a density stratihcation in the inner core. 

Recent [N ii] and [O ill] HST images of NGG 6891 
(Project ID:8390; PI: Hajian) reveal a prolate ellipsoidal 
nebula, as proposed by Guerrero et al. (2000), but neither 
the jet-like structures nor the other kinds of collimated out¬ 
flows are apparent (Eig. 5, upper-left panels). Given that the 
[N ii]/[0 III] line ratio image unveils the faint LISs within the 
much larger scale structures of PNe (Gorradi et al. 1996), we 
constructed the same image using the recent high-resolution 
HST images (Eig. 5, upper-right panel), where no jet-like 
structures are discerned. This provides additional confirma¬ 
tion of our previous conclusion, based on our low-dispersion 
spectra, that there are no jet-like structures in NGG 6891. 
The bottom panels in Eig. 5 display the radial prohle of the 
[N ll] and [O III] emission-line normalized huxes (panels a 
and c) and the [N ll]/[0 ill] line ratio (panels b and d) along 
the slits for the two PAs of 45° and 135°, respectively. 

The ionic and chemical abundances of NGG 6891 are 
presented in Tables 11 and 13. NGG 6891 is a non-type I 
nebula due to its low He and N/0 abundances. The chem¬ 
ical abundances for all nebular components, along both slit 
positions, are found to be the same, within the errors. Our 
chemical abundances of NGG 6891 agree with those derived 
from previous work (Wesson et al. 2005). 

3.5 NGC 6572 

NGG 6572, or PN G034.6+11.8, is a widely studied PN. 
Narrow-band images show an elliptical shell with the ma¬ 
jor axis oriented along PA — 0° and a seemingly toroidal 
structure (Miranda et al. 1999). The [N ll]/[0 III] image ra¬ 
tio highlights the presence of two pairs of knots at PA=15° 
and PA=162° (see Eig. 1 from Miranda et al. 1999). Both 
pairs are located at the tips of two much fainter collimated 
outflows, thus suggesting shock interactions. HST images of 
NGG 6572 do not conhrm the presence of any knotty struc¬ 
ture. Gontrariwise, the HST H(a/[0 III] line ratio image re¬ 
veals a hlamentary structure around the collimated outflows 
indicating strong shock interactions (Guerrero et al. 2013), 
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Figure 5. Upper panels: HST [N ll] image (left-hand panel) and [N ii]/[0 ill] line ratio image (right-hand panel) of NGC 6891. The 
size of the fields shown are 32x32 arcsec^ and 20x20 arcsec^, respectively. The nebular components under analysis are indicated. The 
extraction windows for these components are indicated by the corresponding boxes: the entire NW-SE and NE-SW nebular regions (entire 
NEB: 20.3 arcsec, in red); the NW and SE outer nebular regions (4.9 arcsec, in blue); NE and SW outer nebular regions (6.3 arcsec, in 
cyan); the SE and NW knots (2.1 arcsec, in green); and the inner nebular regions (inner NEB: 4.9 arcsec, in magenta). North is up, and 
East is to the left. Lower panels: The radial profile of the [N ll] and [O ill] emission lines fluxes, normalized to 1.0 and the [N ll]/ [O ill] 
line ratio, along the slits: Panels (a) and (b) correspond to PA=45° whereas (c) and (d) to PA=135°). 


probably responsible for the narrow zones (shock fronts) and 
the enhancement of low-ionization lines. 

The spectral type of NGC 6572 central star is still de¬ 
batable. First, it was classified as a Of-WR type (Mendez 
et al. 1990), then as a WEE star (Tylenda et al. 1993; 
Parthasarathy et al. 1998), and more recently as a [WC]- 
PG 1159 star, by Marcolino et al. (2007). Narrow N III 
AA4634 — 41 and C IV A5808 emission lines (with FWHM of 
15zb2 and 20zb3 A, respectively) are detected in our spectra, 
which is in agreement with the previous WEL-type or [WR] 
classification. 

In order to study the physico-chemical properties of the 
nebula and the pairs of knots, two low-dispersion spectra 
were obtained along the PA=15° and 162°. The emission 
line fluxes as well as c(H/3), Ne and Te are presented in the 
Tables 14 and 16. All the nebular structures used here are 
defined in Fig. 6. The [N ii]/[0 III] ratio profile of NGC 6572, 
along the two slits, are shown in Fig. 6 (panels c and d). 
c(H/3) varies slightly among the nebular component from 


0.24zb0.05 to 0.29zb0.06. These values are lower than those 
published by Liu et al. (2004; 0.48) and Frew et al. (2013; 
0.41), but very close to the value derived by Kwitter et al. 
(2001; 0.32). Regarding the Te, no variations were found be¬ 
tween the inner nebula and the pairs of knots (see Table 14 
and 16). Unlike Te, shows significant variations among 
the components. In particular, the knots are found to have 
lower density than the inner nebular regions, by a factor 
between 2 and 4. Moreover, the high densities derived from 
the Ar and Cl diagnostic line ratios imply a strong density 
stratification in its central region. Ae[S ii] is found to be 
lower than those previously published by Wang et al. (2004) 
and Kwitter et al. (2001), whilst Ae[Ar iv] and Ae[Cl III] are 
higher. The very high of this nebula, close to the upper 
limit, indicate that any small change in the line ratio will 
results in substantially differences values. This can ex¬ 
plain the differences found in among the aforementioned 
works and ours. 

Tables 15 and 17 list the ionic and total abundances 
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Figure 6. Upper panel: HST Hq/[0 hi] image of NGC 6572. The 
size of the field shown is 22x22 arcsec^. The nebular components 
under analysis are indicated in the image. The extraction windows 
for these components are indicated by the corresponding boxes: 
the entire nebular regions (entire NEB: 17.5 arcsec, in red); the 
inner nebular regions (inner NEB: 7.7 arcsec, in cyan); and the 
four, NE, NW, SE and SW knots (2.1 arcsec, in green). North 
is up, and East is to the left. Lower panels: The radial profile of 
the [N ll] and [O ill] emission-line fluxes, normalized to 1.0, and 
the [N ll]/ [O III] line ratio, along the slits. Panels (a) and (b) 
correspond to PA=11°, whereas (c) and (d) to PA=160°). 


for all the nebular components of NGC 6572 and for both 
slit positions. We do not find any significant enrichment of 
N or S, with respect to the rest of the nebula that could 
be associated with the brighter [N ll] and [S ll] lines found 
in the knots. Chemical abundances for a given element are 
the same for all the components. Comparing our chemical 
abundances with those derived by Perinotto et al. (2004) and 
Henry et al. (2004), we turn out that He and N abundances 
agree with those from Perinotto et al. (2004), whilst O, Ne, 
and S abundances better agree with those from Henry et 
al. (2004). There is a substantial difference between these 
two studies. Moreover, the S knot (PA=160°) and the N 
knot (PA=15°) show enriched Ne abundance compared to 
the other two knots. Regarding our analysis, NGC 6572 is a 
non-type I PN, in contrast to the previous classification by 
Henry et al. (2004). 


4 RESULTS: MORPHO-KINEMATIC 
PROPERTIES 

Here we discuss the main morphological and kinematic char¬ 
acteristics of each nebula, and their different components, 
based on previously published works and echelle data taken 
from the SPM KCGPN. The wide range of expansion veloci¬ 
ties of LISs, from a few tens to hundreds km s“^(Gongalves 
et al. 2001), may suggest that shocks play a major role in 
the formation and excitation of these structures. Therefore, 
a parallel study of the morpho-kinematic and the physico¬ 
chemical properties (the emission line ratios, ionic and total 
chemical structure) is necessary to enlighten our knowledge 
of these structures. The de-projected expansion velocities of 
the nebular components, for each PN, are given in Table 1. 

4.1 IC 4846 

A morpho-kinematic analysis of IC 4846, performed by Mi¬ 
randa et al. (2001), reveals the presence of three pairs of 
elongated structures labelled as A1-A2 along PA=54°, and 
B1-B2 / C1-C2 along PA=11°. The authors claim that the 
components A1-A2 correspond to a high collimated precess- 
ing jet. The high-resolution HST Ha image of IC 4846, pre¬ 
sented here for the hrst time, provides further evidence that 
supports this hypothesis (see Fig. 2). 

The projected expansion velocities of NE and SW jets 
range from 36 to 52 km s“^, with an average value of 
44 km s“^ (Miranda et al. 2001). Considering that the incli¬ 
nation angle of the precession axis, with respect to the plane 
of sky, is 13-26°, the de-projected Vexp of the NE and SW 
jets are estimated to be between 100 and 190 km s“^. The 
inner NEB of IC 4846 (or the C1-C2 features in Miranda et 
al. 2001) expands with a velocity of ~20-25 km s“^. 

4.2 Wray 17-1 

The kinematic study of this peculiar nebula has been per¬ 
formed by Corradi et al. (1999). The NW and SE hlaments 
(labelled as A and B in Corradi et al. 1999) show a lin¬ 
ear increase of expansion velocities up to ~30 km s“^ in 
the [N ll] emission line. Incidentally, the NW and SE knots 
at the outer parts of the nebula (labelled as A' and B' re¬ 
spectively in Corradi et al. 1999), show peculiar expansion 
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Table 1. De-projected expansion velocities of LISs and large-scales components. Errors are between 10% and 20% due to the uncertainties 
on the inclination angle. 


PN name 

PA 

(°) 

Jet/Jet-like 
or Filament 
( km s-i) 

Jet/Jet-like 
or Filament 
( km s-i) 

Knot 

( km s-l) 

Knot 

( km s-l) 

Inner NEB 

( km s-l) 

Outer NEB 

( km s-i) 

IC 4846 

54 

^100 (NE) 

^100 (SW) 

— 

— 

20-25 

- 

Wray 17-1 

75 

- 

- 

40-50 (W) 

40-50 (E) 

40 

40“ 

Wray 17-1 

155 

70-80 (NW) 

60-70 (SE) 

>100 (NW) 

>40 (SE) 

40 

- 

K 1-2 

153 

30-35 

- 

30-35 (NW) 

30-35 (SE) 

- 

- 

NGC 6891 

45 

- 

- 

- 

- 

10 

28 

NGC 6891 

135 

- 

- 

80 (NW) 

80 (SE) 

17 

45 

NGC 6572 

15 

- 

- 

46 (NE) 

46 (SW) 

18 

- 

NGC 6572 

169 

- 

- 

65 (NW) 

55-65 (SE) 

18 

- 


® This expansion velocity refers to the E and W arcs. 



s ^hel hm s '] [km s ’] 

Figure 7. Echelle spectra of Wray 17-1 drawn from KCGPN. 
Left panel: [N ll], along PA=155°. Middle panel: [O ill], along 
PA=155°. Right panel: [O iii], along PA=75°. 

velocities compared to the filaments. In particular, the NW 
and SE knots are found to be both red-shifted, with I4xp of 
~60 km and ~10 km s“^, respectively, whilst the NW 
and SE filaments are blue and red-shifted, respectively. On 
the other hand, the E and W knots (respectively C and D, in 
Corradi et al. 1999) show expansion velocities of 15 km s“^ 
with the E knot being blue-shifted and the W knot red- 
shifted. The arc structures that surround these knots show 
similar expansion velocities. 

High-dispersion [O III] PV diagrams along PA=60° and 
135°, were acquired from the SPM KCGPN to better con¬ 
strain the morpho-kinematic structure of Wray 17-1 (Eig. 7). 
These PV diagrams clearly show the presence of an outer 
shell bright in the [O ill] emission, but absent in the [N ll] 
one. This structure expands with a velocity of 60 km s“^, 
with the NW part red-shifted and the SE blue-shifted. The 
LISs seem to interact with the outer shell and indicate shock 
excitation. 


If the central star of Wray 17-1 ejected stellar material 
along a geometrical cone, this could give a good explanation 
for the peculiar velocity field of the filaments and knots. In 
particular, the filaments could be the result of the interac¬ 
tion between the cone and the outer [O III] shell, whilst the 
NW knot, located exactly behind, and being well aligned 
with the filaments. Under this hypothesis, we can explain 
the blue-shifted NW filament and red-shifted NW knots. In 
the case of the SE filament and knot, this scenario implies 
that the two components are related. We calculate the in¬ 
clination of this cone to 20-25° with respect to the line of 
sight. This implies that the de-projected velocities of the fil¬ 
aments and the knots are at least 2.5-3 times higher. On 
other hand, the E and W knots are found to share the same 
expansion velocities of the surrounding nebula, with a de- 
projected velocity ~40-50 km s“^. The inner regions, NEB 
(NE, SW, NW and SE) appear to expand with a velocity of 
~40 km s“^ 


4.3 K 1-2 

K 1-2 shows the presence of a collimated filament and sev¬ 
eral knots oriented along PA=-27° (Corradi et al. 1999). 
These authors pointed out that all component appear to 
move with a velocity similar to that of the main nebula, 
25 km s“^, under the assumption that the NW filament 
(labelled A in Corradi et al. 1999) is located within the main 
spherically symmetric nebula. This assumption implies that 
the filaments and knots are oriented at a minimum incli¬ 
nation angle of 40° with respect to the line of sight. Given 
that the inclination angle of the spectroscopic binary has 
been determined 50° (Exter et al. 2003), the de-projected 
expansion velocity of the knots and jet-like structure would 
be ^30-35 km s“^. 

4.4 NGC 6891 

Guerrero et al. (2000) studied the kinematic properties 
of this nebula and reached the conclusion that a highly- 
collimated outflow (jet) is present along PA=135°. The max¬ 
imum de-projected velocity is 80 km s“^, at the outer edges 
(A-A' component or the pair of knots), considering an in¬ 
clination angle with respect to the plane of sky of 10°. Re- 
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garding the inner and outer NEB, the authors adopted an 
ellipsoidal shell model with a homologous expansion law to 
reproduce the observed data. They concluded that the inner 
NEB expands with a velocity of 10 km s“^ along the minor 
(PA=45°) and 17 km s~^ along the major axis (PA=135°), 
while the outer NEB expands faster with 14xp=28±5 km s“^ 
along the minor (PA 70°) and 45±5 km s“^ along the major 
axis (PA~160°). 

4.5 NGC 6572 

Eollowing Miranda et al. (1999), NGC 6572 shows the pres¬ 
ence of two pairs of low-ionization knots oriented along 
PA=15° and PA=-18°. Moreover, two, much fainter, ellip¬ 
soidal structures can also be discerned in these directions. 
Eig. 6 shows, in fact, that the two pairs of knots are not really 
low-ionization nebular condensations, but resemble those of 
knots due to the low angular resolution of the ground-based 
images. We think that these knots actually result from the 
interaction between the outflows and the outer shell (AGB 
envelope), consistent with theoretical shock models (Dopita 
1997; Raga et al. 2008). 

The projected expansion velocity is 39 km s“^ for the 
NW and the SE knots, and 8 km s“^ for the NE and SW 
knots. It is worth mentioning here that the SE knot shows a 
peculiar decreasing velocity, from 38 to 33 km s“^ (Miranda 
et al. 1999). The NGC 6572 inner NEB shows a projected 
velocity of 14 km s~^. Moreover, two clearly separated max¬ 
ima are discerned in the PV diagrams, indicating the pres¬ 
ence of a ~14 km s“^ toroidal structure at the centre of 
the nebula. This toroidal component has an inclination an¬ 
gle of 38° with respect to the plane of sky (Hora et al. 1990; 
Miranda et al. 1999). In accordance with this inclination an¬ 
gle, the de-projected expansion velocities are calculated to 
18 km s“^ for the torus and the inner NEB, ~65 km s“^ 
for the NW knot, and ~55-65 km for the SE knot. The 
ellipsoidal structure along PA=15° is less tilted, by ~10 °, 
and we obtain a de-projected velocity of ~46 km s“^ for 
both the NE and the SW knots. 


5 DISCUSSION 

5.1 Electron temperatures 

Te was estimated for several individual regions e.g., knots, 
filaments, jets and jet-structures, inner and outer nebular 
regions and arcs, as well as the entire nebulae. We mainly 
used the [O ill], [N ll] and [S ll] diagnostic lines. Te is found 
to range between 9700 K and 14700 K, which is typical of 
photo-ionized nebula, with an average of Te[0 li] of 11500 K 
and Te[N ll] = 11340 K. Moreover, among the various nebular 
components, no significant variation in Te is found, which 
agrees with previous studies (e.g. Gongalves et al. 2009; and 
references therein). 

5.2 Electron densities 

Regarding the electron densities, previous studies on the 
main nebular components (rims and shells) versus LISs have 
shown that the latter have systematically lower densities 
than the former, by a factor between 1 and 10 (see Gongalves 


et al. 2009; and references therein). This result comes into 
conflict with the theoretical predictions, on which the LISs 
are denser than the surrounding medium. 

The most likely explanation for this discrepancy be¬ 
tween models and observations may be that the models are 
referred to the total density of gas (dust, atomic and molecu¬ 
lar as well) while Ne corresponds only to the ionized fraction 
of the gas (Gongalves et al. 2009). In this work, we confirm 
that LISs have the same or even lower densities than their 
surrounding medium. Indeed, in a paper in prep, we report 
our first successful detection of the molecular hydrogen emis¬ 
sion, at 2.12 /xm, associated to the LISs of two PNe (K 4-47 
and NGG 7662). Although, this detection was tried before 
by the group, the fact that LISs have in general smaller 
scales than the main bodies did not help. It seems that only 
8-m class telescopes can successfully detect the H 2 associ¬ 
ated with LIS in Galactic PNe. See Manchado et al. (2015) 
for a similar work, though not associated with LISs directly, 
but with the molecular torus of NGG 2346. Moreover, it 
turns out that the molecular torus of this PN is composed 
of knots and filaments. On the other hand, in the Helix neb¬ 
ula (Meaburn et al. 1992; ODell & Handron 1996; Huggins 
et al. 2002), which posses the best-known LISs of PNe, the 
association between LISs and the molecular gas is known for 
years. 

Images of PNe in molecular hydrogen (H 2 ) and optical 
emission lines (e.g. [N ll], [S il] and [O l]) display similar 
morphologies (Guerrero et al. 2000; Bohigas 2001). In ac¬ 
cordance with the theoretical work by Aleman & Gruenwald 
(2011), the peak intensities of these optical lines and the H 2 
v=l-0 co-rotational lines occurs in a narrow transition zone 
between the ionized and the neutral (photo-dissociation) re¬ 
gions. It is known that molecular hydrogen in PNe can be 
excited either by shocks or by the absorption of UV-photons 
(fluorescence) emitted by the central star. Therefore, the de¬ 
tection of strong [N il] and [O l] emission lines in LISs may 
be the result of shock excitation, due to the moderate ex¬ 
pansion velocities (see § 4). 

Shock models can successfully reproduce the strong low- 
ionization lines from LISs (K 4-47; Gongalves et al. 2004, 
Raga et al. 2008). In particular, Raga et al. (2008) reached 
the conclusion that the lower the photo-ionization rate, the 
higher the Ne diagnostic line ratios, which then resemble 
the typical values of shock-excited regions. Besides the ex¬ 
pansion velocities of LISs, the evolutionary age and photo¬ 
ionization rate of the central source are also key parame¬ 
ters for probing the excitation and formation mechanisms of 
LISs, thus unveiling their true nature. 

5.3 Ionic and total abundances 

The chemical abundances were calculated for the several 
nebular components used here for the analysis of the 5 
PNe. Gomparison with those abundances derived in previous 
studies shows a reasonable agreement with previous results, 
not only for the entire nebulae but also for each component 
individually. 

Eor a given PN, and element, all the components ex¬ 
hibit the same chemical abundances, within the uncertain¬ 
ties. This implies that the strong low-ionization lines, such 
as [N ii], [S li] or [O l] of LISs, may not be the result of 
an overabundance of these elements in the nebula. Eig. 8 
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Figure 8. Abundance patterns: log(N/0) vs. He/H (left-hand panel) and log(N/0) vs. log(N/H)+12 (right-hand panel). Red pentagons 
correspond to LISs, green pentagons to jets, filaments and jet-like structures. Nebular regions (inner and outer NEBs) are cyan filled 
diamonds whilst the entire NEBs are black filled diamonds. The dashed line represents the best linear fit to our data (see text). 


displays log(N/0) vs. log(N/H) + 12 (left-hand panel) and 
log(N/0) vs. He/H (right-hand panel) for all the PNe in 
our sample and their components. A linear relation be¬ 
tween log(N/0) and log(N/H) + 12 is found, which is con¬ 
sistent with the general picture for PNe (Perinotto et al. 
2004). The best linear fit to our data yields the relation, 
log(N/O)=0.74*(log(N/H) + 12)-6.5, with (the goodness- 
of-fit) equals to 0.88, and it is in very good agreement with 
the relation derived from a sample of PNe with WR and 
WEL-type central stars (Garcia-Rojas et al. 2013). Only E 
knot in Wray 17-1 appears as an outlier (easily identifiable 
in the figure), so not considered in the fitting process. The 
log(N/0) vs. He/H plot does not show any correlation be¬ 
tween the LISs and the other nebular components, whereas 
both show He/H<0.14 and log(N/O)<-0.3. These figures, 
following Kingsburgh & Barlow (1994) criteria, implies that 
all the PNe in our sample are non Type-I nebulae. 

In Eig. 9 we plot log(X/0) vs. log(X/H)+12 for X=Ne, 
S and Ar, and clear linear trends are found. We have to 
mention that the apparent distinct regions in these dia¬ 
grams do not really reflect chemical abundance differences 
between LISs and other nebular components (NEBs), but 
potential differences between the PNe in our sample and 
their central stars. Eor instance the LISs with Ne/0>- 
0.55 and log(Ne/H) + 12>8.1 (as well as Ar/0>-2.5 and 
log(Ar) + 12>6.1) correspond to the LISs in K 1-2 and Wray 
17-1, for which we were able to estimate their Ne and Ar 
abundances, whereas the NEBs are those from IC 4846, 
NGC 6572 and NGG 6891. Thus, K 1-2 and Wray 17-1 are 
likely richer, in Ne and Ar, than the rest of the PNe in our 
sample. Surprisingly, LISs in the former PNe are highly- 
ionized, with He being mainly doubly ionized (see Tables 5, 
7 and 9). As for S, all the components (LISs and NEBs) are 
well mixed and cover the same range of values. Trends we 
found in Eigure 9 are those usually seen in PNe between the 
a-element to oxygen ratios. 

In conclusion, the similarities in the chemical abun¬ 
dances among the components in each PN strongly sug¬ 
gest that the intense low-ionization emission lines are the 
result of their excitation mechanisms rather than the chem¬ 
ical overabundance of any element in the LISs. 

Moreover, though a deep discussion about the abun¬ 
dance discrepancy factor (adf) is out of the scope of the 
present paper, it is relevant to note that adjs, close binaries 
as PN central stars (so, post GE stars) and the presence of 


LIS could somehow be related. As briefly discussed in the In¬ 
troduction, the formation LIS, and particularly that of jets 
older than the main shells, is expected in the case of post 
common-envelope PN central stars (Gongalves et al. 2001; 
Miszalski et al. 2011). Very recently, it was shown that the 
PNe having the highest adjs (10 to 300) all have post GE 
central stars (Gorradi et al. 2015). In addition of confirm¬ 
ing the mixing of the bi-chemistry in these PNe (hot normal 
component plus hydrogen-deficient cold inclusions; Liu et al. 
2004), Gorradi et al. (2015) point out the present difhcult 
in depicting an evolutionary scenario consistent with all the 
observed properties of the highest ad/PNe. Some of the sce¬ 
narios these authors discuss could indeed be related with 
the presence of LISs via, for instance, the photo-ionization 
of the neutral material that was in the orbital plane during 
the GE phase. 


5.4 Excitation mechanism 

Besides the emission line ratios, kinematic information -in 
fact expansion velocities- is essential to fathom the forma¬ 
tion and excitation mechanisms of LISs in PNe. The diag¬ 
nostic diagrams, based on the shock models from Raga et al. 
(2008), are used here in order to perform a more comprehen¬ 
sive study on the LISs and the other nebular components by 
combining their emission lines and expansion velocities. 

In particular, Raga et al. (2008) modelled the evolu¬ 
tion of a high-speed (I00-I50 km s“^) and high-density 
(10^ cm“^) knot moving away from a photo-ionizing source 
through a low-density medium, for a range of initial condi¬ 
tions. The simulated knot is a spherical cloudlet with tem¬ 
perature of 10000 K, and size of 10® cm. The surround¬ 
ing medium was assumed to be uniformly filled with a low- 
density gas of 10^ cm“^ and the same temperature (10000 K) 
of the knot. A black-body approximation was assumed for 
the central photo-ionizing source, with a L=5000 L0 and 
two Teff of 50000 K and 70000 K, which correspond to an ion¬ 
izing photon rate of iS^^d.dlx 10^^ s“^ and 4.66x10^^ s“^, 
respectively. Moreover, the effect of the distance between 
the knot and the central source was also studied for three 
specific distances of 3, I, 0.3x10^® cm. 

The emission line ratios of the knots, for two differ¬ 
ent epochs of the simulations (200 yr in red and 400 yr, 
in green), are presented in Eig. 10, superposed to the ob- 
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Figure 10. Diagnostic diagrams distinguishing the shock-excited and the photo-ionized regions. Panels a, b, c, and d are the diagnostic 
diagrams from Raga et al. (2008), on which we superimpose the results for our sample of PNe. The regions with diagonal solid lines 
represents the loci of fast LISs, whereas the regions with dashed lines represent the location of photo-ionized structures (e.g. inner and 
outer NEBs that are in fact shells and rims). Panels e and f display the diagnostic diagrams by Sabbadin et al. (1977), with the density 
ellipses introduced by Riesgo & Lopez (2006). 


served emission line ratios from this work. Unfortunately, 
these models do not provide any predictions for the prop¬ 
erties of the other nebular components like shells, rims 
and haloes. The symbols in Fig. 10 correspond to differ¬ 
ent expansion velocities, Tgff, and distances to the central 
source, -boxes: Vexp=100 km s“^ and Teff=50000 K (mod¬ 
els a/b/clOO); crosses: Ueccp=100 km s“^ and Tgff=70000 K 
(models A/B/ClOO); and x symbols: Veccp=150 km s“^ 
and Teff=70000 K (model A/B/C150). Letters from “a” to 
“c” (and “A” to “C”) indicate the different distances from 
the central source and thus the different photo-ionization 
rates (small symbol = small distance = closer to the cen¬ 
tral source). The solid and dash-dotted lines enclose the 
regime of fast LISs (e.g. knots, jets) and rims and shells (or, 
more general, photo-ionized regions) taken from Raga et al. 
(2008). The open pentagons correspond to the LISs (knots: 
red and jets/jet-likes/filaments: green), the filled cyan dia¬ 
monds correspond to the inner and outer NEB regions, and 
the filled black diamonds to the entire NEB. Moreover, the 


larger the points’ size the higher the de-projected expansion 
velocity (see Table 1). The velocities are normalized to the 
average expansion velocity of PNe, 31 km s“^ (Pereyra et 
al. 2013). Eor the entire NEB (black symbols), we consid¬ 
ered an expansion velocity equal to the value just quoted, 
so emphasizing the different expansions among the nebular 
components in the diagnostic diagrams. No trend is found 
between the emission line ratios and Vexp of the LISs and 
the nebular regions. This suggests that any attempt to disen¬ 
tangle the photo-ionized and shock-excited regions, or even 
investigate their excitation mechanisms, requires the evolu¬ 
tionary age and photo-ionization rates of the central star, 
apart from the expansion velocities themselves. 

All LISs in this work are found to be located in the 
regime of fast LISs, except from the pair of knots in 
NGC 6891. The latter knots are displaced from the region of 
fast LISs, in each diagram, despite their high expansion ve¬ 
locities. Moreover, the NE and SW outer NEBs in NGC 6891 
(the jet components in Guerrero et al. 2000) also lie in the 
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Figure 11. Plots of [N ii]/Hq; (panel a), [O Ilj/Ho; (panel b), [O ll]/[0 ill] (panel c) and [S llj/Ha (panel d) vs. log(fs/ioc/es/f*)- The 
vertical strips represent log(fs/j,ocA:s/f*) = 0. K 1-2 and Wray 17-1 structures were arbitrarily placed into the strips, since there are no 
distance measurements for these PNe. 


regime of photo-ionized regions and close to the inner NEBs 
of the nebula (NE and SW outer NEBs). This is a further 
indication that there is no difference among the regions ana¬ 
lyzed in NGC 6891. We already pointed out that there are no 
jets in NGC 6891 (see §3.4). Knots, on the other hand, show 
low-ionization lines slightly enhanced. We believe that this 
may be associated with the low contribution of the photo¬ 
ionization in these knots (higher distance from the central 
star), resulting in more efficient shock excitation. 

Adopting a distance of 2.5 kpc for NGC 6891 (Phillips 
2005) and an angular distance from the central star of 
9.1 arcsec, the linear distance is ~3.4xl0^^ cm. Given that 
the Teff and L of the central star are 50000 K (Mendez 
et al. 1988) and 10000 L© (Mendez et al. 1992; re-scaled 
to the distance used here), respectively, the ionizing pho¬ 
ton hux (E) reaching the knots is fjH=5.46x 10^^ cm“^ 
s“^. NGC 6891’s knots he closer to the Raga’s ClOO model 
(Eeccp=100 km s“^, E)=3xl0^^ cm and E=3.73x10^1 cm“^ 
s“^) than to the clOO model (Wccp^lOO km s“^, 
D=3xl0^^ cm and E=2.73x10^1 cm“^ s“^). The main dif¬ 
ference between these two models is the ionizing photon hux, 
which depends on the Tgff. Models “C” are those with the 
smallest distance from the central stars and therefore have 
higher ionization rates. Hence, NGC 6891 appears to be 
a nebula with a higher ionization rate (so, mainly photo- 
ionized) whilst the contribution of shocks must be very low. 
Nevertheless, the high expansion velocities of the pair of 


^ The ionizing photon hux, E, is calculated by means the equa¬ 
tion fi,=Si,4:7r d?^ where S'* is the photon rate and d is the distance 
of the cloudlet from the central star. 


knots (~80 km s ^), results in the enhancement of the low- 
ionization lines. 

As for NGC 6572, the inner NEBs are in a transition 
zone between the rims/shehs and fast LISs regimes, whilst 
the two pairs of knots are well placed on the regime of the 
fast LISs. The de-projected expansion velocities of the two 
pair of knots are at least 2.5-3 times higher than the ex¬ 
pansion velocities of the inner NEBs, then suggesting higher 
contribution from shocks. The angular distances of the knots 
from the central star are 6.67 arcsec (NE and SW), 7.5 arcsec 
(NW) and 8.4 arcsec (SE). Considering a distance of 1.4 kpc 
(Mehema 2004), we obtain linear distances from the central 
star of 1.35, 1.5 and 1.71 x 10^^ cm. These distances give 
ionizing photon huxes between 2.3 and 3.4x10^° cm“^ s“^, 
considering Tgff=60000 K and L=2250 L© (Shaw & Kaler 
1989; re-scaled to distance adopted in this work). This ioniz¬ 
ing photon huxes imply a low contribution in the excitation 
of the knots due to the absorption of UV photos. Therefore, 
the enhancement of the low-ionization lines in these struc¬ 
tures should be a result of the shocks. The higher ionizing 
photon hux of the inner NEBs, due to their smaller distance 
to the central star, in conjunction with their much lower ex¬ 
pansion velocities, explain their position on the diagnostic 
diagrams in the regime of photo-ionized regions. 

Unlike the knots of NGC 6572, the jets of IC 4846 are 
located in the transition zone, despite their high expansion 
velocities (>100 km s“^). Eor Tgff=50000 K and L=2067 L© 
(Shaw & Kaler 1989; rescaled to the distance used here), 
the ionizing photon rate is 1.36 x 10^^ s“^, and for an an¬ 
gular distance from the central of 3.5 arcsec and D—3.72 
kpc (Phillips 2004), the linear distance is then 3x10^^ cm. 
This implies an f,=1.20x10^^ cm“^ s“^, which is approxi¬ 
mately three times lower than that of ClOO model. Despite 
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log(S/H) + 12 



Figure 9. Abundance patterns: log(Ne/0) vs. log(Ne/H)+12 
(upper panel); log(S/0) vs. log(S/H)+12 (middle panel) and 
log(Ar/0) vs. log(Ar/H)+12 (lower panel). Symbols are the same 
as in Fig. 8. A trend of increasing Ne/O, S/O and Ar/O with in¬ 
creasing Ne, S and Ar, respectively, is found. 


this low value, IC 4846 exhibits stronger [N il], [S ii] and 
[O l] emission lines than in model Cl00, which means that 
shock interactions are responsible for this enhancement. 

All LISs of Wray 17-1 are well placed in the regime of 
fast LISs, with intense low-ionization lines. Unfortunately, 
there is no distance measurement for this object. Neverthe¬ 
less, assuming a linear distance of 3x10^^ cm between the 
knots and the central star (the smallest distance among the 
shock models), and an angular distance of 43 arcsec from 
the central source, a D =0.46 kpc is obtained, which seems 
unrealistic. This may suggest that the knots are at a larger 
distance from the central star, e.g. 3x10^® cm (or D ~5 
kpc). We believe that the strong low-ionization lines found 
in this nebula are resulting from shocks, whilst the contri¬ 
bution of the stellar UV photons may be very low. The ob¬ 
served spectra from the inner NEBs are in agreement with 
this argument. 

Regarding the H(a/[N il] vs. Hq:/[S ii] and [S ll] 
AA6716/6731 vs. Hq;/[S ii] diagnostic diagrams, introduced 


by Sabbadin et al. (1977), it is clear that all shock-excited 
features are located at the bottom-left region whilst the 
photo-ionized features at the upper-right region (lower pan¬ 
els in Fig. 10). The zone between 0.2< log(H(a/[N ll])<0.8 
defines the transition zone where both excitation mecha¬ 
nisms are efficient. In general, hot and dim central stars are 
able to reproduce emission lines ratios that resemble those 
of shock-excited nebulae without necessarily being shocked- 
excited (Raga et al. 2008). This could explain the pres¬ 
ence of some elliptical and round PNe in the region with 
log(HQ;/[N ii])< 0 and log(HQ;/[S ii]])< 0.4 (Sabbadin et. 
1977; Riesgo & Lopez 2006). 

In accordance to the previous analysis, besides the emis¬ 
sion line ratios, kinematic information is required in order 
to probe the excitation mechanisms of the LISs in PNe. 
Hence, a more general approach, which would take into ac¬ 
count the contribution of shock- and photo-ionization is ac¬ 
tually required here. We quantitatively investigate this by 
defining the dimensionless parameter log(fs/iocfcs/f+), where 
fshocks is the ionizing photon flux in units of cm“^ s~^ for a 
shock with W, and U is the stellar ionizing photon flux at 
a given distance, fshocks is as follows: we calculate the total 
energy flux emitted due to the shock, by using the formula 
Fshocks=2.28x 10“^(U/100 km s“^)^ x(ne/cm“^) in units 
of erg cm“^ from Dopita & Sutherland (1996). Then, 
we convert this energy flux (Fshocks) to photon flux (fshocks), 
by dividing by the average energy of a photon (L/S+) emit¬ 
ted from the central star. This ratio should give us of the 
contribution of each mechanism to the excitation of the gas. 

In Fig. 11 we present the new diagnostic diagrams 
[N ii]/Ha, [S ii]/Ha, [O ii]/Ha and [O ii]/[0 ill] vs. 
\og{fshocks/U) that allows us to distinguish nebular com¬ 
ponents with different excitation mechanisms. It is notice¬ 
able that there are two distinct regimes in these diagrams: i) 
the upper-right corner with the shock-excited structures and 
ii) the lower-left corner with the photo-ionized structures. 
A transition zone, where both excitation mechanisms are 
important, can also be considered. All the nebular regions 
(NEB) lie in the lower-left region, with \og{ishocks/U)<-‘^, 
whilst the knots and jets lie in the upper-right part of the 
diagrams, with log(fs/iocA:s/f+)>-l- Only the pair of knots 
in NGC 6891 show \og{ishocks/U)<-'^, consistently with our 
previous conclusion that they must be mainly photo-ionized 
with a low contribution from shocks. The models a/blOO, 
A/BIOO and A/B150 also lie in the shock-excited regime 
(upper-right corner), whilst the models clOO, ClOO and C150 
are found in the transition zone. The latter models have 
the higher ionizing photon flux among all the models, which 
means that the excitation due to the absorption of stellar UV 
photons is more efficient. Any decrease of Wxp, for instance 
from ClOO to C70 or C40, will result in further decrease of 
\og{fshocks/U) and will eventually place these models closer 
to the regime of photo-ionized regions. 

The points representing K 1-2 and Wray 17-1 have been 
arbitrarily added to the diagram, at log(fs/iocA:s/f+)=0, since 
there are no distance measurements for these objects and 
the U cannot be calculated. Based on their emission line ra¬ 
tios, one can assume that the jets and knots of these two 
objects have, indeed, log{ishocks/U) close to 0. Assuming 
log(fs/ioc/es/U)=0 for the knots and jets, log(fs/iocfcs/U)=-l 
for the arcs, and L=5000 LO, we obtain a rough estimation 
of their distances. Wray 17-1 should be located between 3 


© 0000 RAS, MNRAS 000, [^T? 





















































16 Akras & Gongalves 


and 7 kpc, where for K 1-2, we obtain a very high distance, 
>20 kpc. This suggests that the lumonosity of the central 
star in K 1-2 is lower than 5000 L© and/or the expansion 
velocities of the knots are higher than 35 km s“^. For in¬ 
stance, a L=2000 L© and Veccp=70 km s“^ gives us a more 
reasonable distance of ~5 kpc. 


5.5 Central Star 

Despite the various attempts to explain the formation of col¬ 
limated (those that appear in pairs) LISs, their origin still 
remains debatable (see in Gongalves et al. 2001). The most 
prominent mechanisms is the interaction of a close binary 
system due to the mass-transfer exchange. The formation 
of an circumstellar/circumbinary disk via binary interaction 
provide those conditions to form bipolar outflows and/or 
collimated jets (Blackman et al. 2001, Nordhaus & Black¬ 
man 2006). Some models suggest that the formation of jets 
occurs after the common envelope (CE) phase of a close bi¬ 
nary star (Soker & Livio 1994). Pairs of jets with kinematical 
ages much younger than the main nebular structures have 
indeed been found in some post-CE PNe, like NGC 6778 and 
NGG 6337 (Tocknell et al. 2014). Although, the majority of 
the post-GE PNe are found to exhibit jets that predate the 
main nebular structures, suggesting that the polar ejection 
occurs before the GE phase (see Tocknell et al. 2014; Jones 
et al. 2014). The discovery of several new close binary stars 
in PNe has provided strong evidence of a link between the 
close binary post-common envelope PNe and LISs (Miszalski 
et al. 2009b). 

Moreover, a possible link between the spectral type of 
the PN nucleus and LISs has also been pointed out by Mis¬ 
zalski et al. (2009a). According to these authors, more than 
half of PNe that are known to possess LISs (~65%) have an 
emission line type central star (e.g. Of, [WR] or PG1159). In¬ 
cidentally, all the PNe studied in this work have an emission¬ 
line star, except K 1-2. In particular, two of them have a 
WEL type nucleus (IG 4846, NGG 6891), one has a [WG]- 
PG1159 (NGG 6572), and one a PG-1159 nucleus (Wray 17- 
1). The central star of K 1-2 has been classified as a K2 V 
but it is a known post-GE close binary system (Exter et al. 
2003; De Marco 2009). 


6 CONCLUSIONS 

The properties like Ae, Te, ionic and total chemical abun¬ 
dances were calculated for various nebular components, like 
nebular regions (inner and outer NEBs) and LISs (jets, jet¬ 
like, filaments and knots), as well as for the entire nebu¬ 
lae of a small sample of PNe (IG 4846, Wray 17-1, K 1-2, 
NGG 6891 and NGG 6572). We confirmed that LISs are 
less dense (in terms of Ne) than the surrounding medium, 
whereas their Tg are the same as that of the other compo¬ 
nents of a given PNe. 

There is no chemical abundance enrichment in LISs that 
could be attributed to the strong low-ionization lines. Eor a 
given PN, all the nebular components have the same chemi¬ 
cal abundances within the errors, reflecting on the chemical 
composition of the entire nebula as well as the evolution of 
the central star. 


A comparison of between the emission lines of the nebu¬ 
lar components and those predicted by the the shock models 
reveals that LISs are predominantly shock-excited regions. 
This shock-excitation of LISs results in the enhancement 
of the low-ionization lines, whereas the other nebular com¬ 
ponents are mainly photo-ionized, so less prominent in the 
low-ionization emission lines. 

In order to probe the contribution of shocks and UV 
photons in the excitation of the LISs, we dehned the param¬ 
eter log(fs/ioc/es/f*), where ishocks and f* are the ionizing pho¬ 
ton fluxes due to the shock interactions and the central star 
continuum, respectively. We then proposed new diagnostic 
diagram that involves the emission-line ratios [N il]/H(a, 
[O li]/Ha, [S ii]/Hq; and [O ii]/[0 ill] versus the new excita¬ 
tion parameter, log(fs/iocfcs/G)- A robust relation differenti¬ 
ating the regimes of shock-exited (log(fs/iocA:s/R)>-l) from 
that of the photo-ionized components (log(fs/iocfcs/R)<-2), 
in addition to a transition zone (-2< log(fs/iocfcs/R)<-l), 
where both mechanisms are equally important, is found. 
Most of the LISs analyzed in this work were found to be well 
located in the regime of shock-excited structures, apart from 
the pairs of jets and knots in NGG 6891. The former was 
found to be photo-ionized as much as its surrounding neb¬ 
ula, while the latter is placed in the transition zone, which 
indicates a contribution from shocks. All the inner and outer 
NEBs are well placed in the region of the photo-ionized com¬ 
ponents. 


ACKNOWLEDGMENTS 

We thank, enormously, R. Gorradi and H. E. Schwarz (in 
memoriam) who provide us the Danish telescope data we use 
in the present paper. We also would like to thank the anony¬ 
mous referee for his/her detailed report on the manuscript, 
which indeed helped us to significantly improve the paper. 
The work of SA is supported by GAPES post-doctoral fel¬ 
lowship ‘Young Talents Attraction’ - Science Without Bor¬ 
ders, A035/2013. This work was also partially supported 
by EAPERJ’s grant E-26/111.817/2012 and GAPES’s grant 
A035/2013. The analysis here presented is based on obser¬ 
vations obtained at the 2.5-m INT (IDS), operated by the 
Isaac Newton Group on the island of La Palma in the Span¬ 
ish Observatorio del Roque de los Muchachos of the Institute 
de Astrof^ sica de Ganarias. Our work is also based on the 
1.54-m Danish telescope (DEOSG) data, operated by the 
European Southern Observatory (ESO) at La Silla, Ghile, 
and on the 2.1-m SPM telescope, operated by the UNAM at 
the Observatorio Astronomico at San Pedro Martir Obser¬ 
vatory in Baja Galifornia, Mexico. We also use NASA/ESA 
HST data, obtained at the Space Telescope Science Insti¬ 
tute. 


REFERENCES 

Acker A., Raytehev B., Stenholm B., Tylenda R., 1991, 
A&AS, 90, 89 

Acker A., Marcout J., Ochsenbein E., Stenholm B., 
Tylenda R., Schohn G., 1992, The Strasbourg-ESO Gata- 
logue of Galactic Planetary Nebulae. Parts I, 11. 

Aleman L, Gruenwald R., 2011, A&A, 528A, 74 


© 0000 RAS, MNRAS 000, [^T? 


Low-ionization structures in planetary nebulae 17 


Balick B., Preston H. L., Icke V., 1987, AJ, 94, 1641 
Balick B., Rugers M., Terzian Y., Chengalur J. N. 1993, 
ApJ, 411, 778 

Balick B., Perinotto M., Maccioni A., Terzian Y., Hajian 
A. R., 1994, ApJ, 424, 800 

Balick B., Alexander J., Hajian A. R., Terzian Y., Perinotto 
M., Patriarchi P., 1998, AJ, 116, 360 
Balick B., Frank A., 2002, ARA&A, 40, 439 
Blackman E. G., Frank A., Welch C., 2001, ApJ, 546, 288 
Boffin H. M. J., Miszalski B., Rauch T., Jones D., Corradi 
R. L. M., Napiwotski R., Day-Jones A. C., K oppen J., 
2012, Science, 338, 773. 

Bohigas J., 2001, RMxAA, 37, 237 

Boumis P., Paleologou E. V., Mavromatakis F., Papamas- 
torakis J., 2003, MNRAS, 339, 735 
Boumis P., Akras S., Xilouris E. M., Mavromatakis F., Ka- 
pakos E., Papamastorakis J., Goudis G. D., 2006, MN¬ 
RAS, 367, 1551 

Gahn J. H., Kaler J. B., Stanghellini L., 1992, A&AS, 94, 
399 

Gorradi R. L. M., Manso R., Mampaso A., Schwarz H. E., 
1996, A&A, 313, 913 

Gorradi R. L. M., Perinotto M., Villaver E., Mampaso A., 
Gon calves D., 1999, AJ, 523, 721 
Gorradi R. L. M. Garcia-Rojas J., Jones D., Rodriguez-Gil 
R, 2015, ApJ, 807, 99 

Delgado-Inglada G., Rodriguez M., Peimbert M., 
Stasihska, G., Morisset G., 2015, MNRAS, 449, 1797 
Dopita M. A., Sutherland R. S., 1996, ApJS, 102, 161 
Dopita M. A., 1997, ApJ, 485, L41 
De Marco O., 2009, PASP, 121, 316 
Exter K. M., Pollacco D. L., Bell S. A., 2003, MNRAS, 341, 
1349 

Fitzpatrick E.L., 1999, PASP, 111, 63 
Frew D. J., Bojicic 1. S., Parker Q. A., 2013, MNRAS, 431, 
2 

Garcia-Rojas J., Pena M., Morisset G., Delgado-Inglada 
G., Mesa-Delgado A., Ruiz T. M., 2013, 2013, A&A, 
558A, 122 

Garcia-Segura, G., Langer, N., Rozyczka, M., Franco J., 

1999, ApJ, 517, 767 

Garcia-Segura G., Villaver E., Langer N., Yoon S.-G., Man- 
chado A., 2014, AJ, 783, 74 

Gongalves D. R., Gorradi R. L. M., Mampaso A., 2001, 
ApJ, 547, 302 

Gongalves D. R., Gorradi R. L. M., Mampaso A., Perinotto 
M., 2003, ApJ, 597, 975 

Gongalves D. R., 2004, ASP Gonference Proceedings, 313, 
216 

Gongalves D. R., Mampaso A., Gorradi R. L. M., Perinotto 
M., Riera A., Lopez-Martm, 2004, MNRAS, 355, 37 
Gongalves D. R., Ercolano B., Garnero A., Mampaso A., 
Gorradi R. L. M., 2006, MNRAS, 365, 1039 
Gongalves D. R., Mampaso A., Gorradi R. L. M., Quireza 

G. , 2009, MNRAS, 398, 2166-2176 

Gorny S. K., Schwarz H. E., Gorradi R. L. M., van Winckel 

H. 1999, AVAS, 136, 145 

Gorny S. K., Ghiappini G., Stasihska G., Guisinier F., 2009, 
AVA, 500, 1089 

Guerrero M. A., Miranda L. F., Manchado A., Vazquez R., 

2000, MNRAS, 313, 1 

Guerrero M. A., Toala J. A., Medina J. J., Luridiana V., 


Miranda L. F., Riera A., Velazquez P. F., 2013, AVA, 
557A, 121 

Guerrero M. A., De Marco O., 2013, AVA, 553A, 126 
Jones D., Santander-Garcia M., Boffin H. M. J., Miszalski 
B., Gorradi R L. M., 2014, in Asymmetrical Planetary 
Nebulae VI Gonference, 43. 

Henry R. B. G., Kwitter K. B., Balick B., 2004, AJ, 127, 
2284 

Hora J. L., Hoffmann W. F., Deutsch L. K., Fazio G. G., 
1990, ApJ, 353, 549 

Kingsburgh R. L., Barlow M. J., 1994, MNRAS, 271, 257 
Kaler J. B., 1970, ApJ, 160, 887 
Kwitter K. B., Henry R. B. G., 2001, ApJ, 562, 804 
Kwok, S., Purton, G. R., Fitzgerald, P. M. 1978, ApJ, 219, 
L125 

Leal-Ferreira M. L., Gongalves D. R., Monteiro H., 
Richards J. W., 2011, MNRAS, 411, 1395 
Liu Y., Liu X.-W., Luo S.-G., Barlow M. J., 2004, MNRAS, 
353, 1231 

Liu X.-W., Storey P. J., Barlow M. J., Danziger 1. J., Gohen 
M., Bryce M., 2000, MNRAS, 312, 585. 

Lopez J. A., Vazquez R., Rodriguez L. F., 1995, ApJ, 455, 
L63 

Lopez J. A., Richer M. G., Garcia-Diaz Ma. T., Glark D. 
M., Meaburn J., Riesgo H., Steffen W., Lloyd, M., 2012, 
RevMexAA, 48, 3 

Manchado A., Stanghellini L., Villaver E., Garcia-Segura 
G., Shaw R. A., Garcia-Hernandez D. A., 2015, ApJ, 808, 
115 

Manchado, A., Guerrero, M. A., Stanghellini, L., Serra- 
Ricart, M. 1996, The I AG Morphological Gatalog of 
Northern Galactic PNs (La Laguna : lAG) 

Marcolino W. L. F., de Araujo F. X., 2003, AJ, 126, 887 
Marcolino W. L. F., de Araujo F. X., Junior H. B. M., 
Duarte E. S., 2007, AJ, 134, 1380 
McGarthy J. K.m Mould J. R., Mendez R. H., Kudritzki 
R. P., Husfeld D., Herrero A., Groth H. G., 1990, ApJ, 
351, 230 

Meaburn J., Walsh J. R., Glegg R. E. S., Walton N. A., 
Taylor D., Berry D. S., 1992, MNRAS, 255, 177 
Mellema G., Eulderink F., Icke V. 1991, AVA, 252, 718 
Mellema G., 1995, MNRAS, 277, 173 
Mellema G., Raga A. G., Ganto J., Lundqvist P., Balick 
B., Steffen W., Noriega-Grespo, A., 1998, AVA, 331, 335 
Mellema G., 2004, AVA, 416, 623 

Mendez R. H., Herrero A., Manchado A., 1990, AVA, 229, 
152 

Mendez R. H., Kudritzki R. P., Herrero A., 1992, AVA, 
260, 329 

Mendez R. H., Kudritzki R. P., Herrero, A., Husfeld D., 
Groth H. G., 1988, AVA, 190, 113 
Miranda L. F., Vazquez R., Gorradi R L. M., Guerrero M., 
A.; Lopez J. A., Torrelles J. M., 1999, ApJ, 520, 714 
Miranda L. F., Guerrero M. A., Torrelles J. M., 2001, MN¬ 
RAS, 322, 195 

Miszalski B., Acker A., Parker Q. A., Moffat A. F. J, 2009a, 
AVA, 505, 249 

Miszalski B., Acker A., Moffat A. F. J., Parker Q. A., Udal- 
ski A., 2009b, AVA, 496, 813 
Miszalski B., Jone, D., Rodriguez-Gil P., Boffin H. M. J., 
Gorradi R. L. M., Santander-Garcia M. 2011, AVA, 531, 
158 


© 0000 RAS, MNRAS 000, [^T? 


18 Akras & Gongalves 


Monteiro H., Gongalves D. R., Leal-Ferreira M. L., Corradi 
R. L. M., 2013, A&A, 560, 102 
Nordhaus J., Blackman, E. G., 2006, MNRAS, 370, 2004 
ODell G. R., Handron K. D., 1996, AJ, 111, 1630 
Palen S., Balick B., Hajian A. R., Terzian Y., Bond H., E., 
Panagia N., 2002, AJ, 123, 2666 
Parker Q. A., Acker A., Erew, D. J., et ak, 2006, MNRAS, 
373, 79 

Parthasarathy M., Acker A., Stenholm B, 1998, A&A, 
329L, 9 

Pereyra M., Richer M. G., Lopez J. A., 2013, ApJ, 771, 114 
Perinotto M., 2000, Ap&SS, 274, 205 
Perinotto M., Morbidelli L., Scatarzi A., 2004, MNRAS, 
349, 793 

Phillips, J. R, 2004, MNRAS, 353, 589 
Phillips, J. R, 2005, MNRAS, 362, 847 
Raga A. G., Riera A., Mellema G., Esquivel A., Velazquez 
R E., 2008, A&A, 489, 1141 

Rauch T., Werner K. 1997, in The Third Gonference on 
Eaint Blue Stars, ed. A. G. D. Philip, J. Liebert, R. A. 
Saffer (Schenectady, NY: L.Davis Press), 217 
Riesgo H., Lopez J. A., 2006, RMxAA, 42, 47 
Sabbadin E., Minello S., Bianchini A. 1977, AVA, 60, 147 
Sabin L., Parker Q. A., Gorradi R. L. M., et ah 2014, 2014, 
MNRAS, 443, 3388 

Sahai R., Morris M. R., Villar G. G., 2011, AJ, 141, 134 
Schwarz H. E., Gorradi R. L. M., Stanghellini L. 1993, in 
lAU Symp. 155, Planetary Nebulae, ed. R. Weinberger & 
A. Acker (Dordrecht: Kluwer), 214 
Shaw R. A., Kaler J. B., 1989, ApJS, 69, 495 
Shaw R. A., Dufour R. J., 1995, PASP, 107, 896 
Shaw R. A. 2012, in lAU Symp. 283, Planetary Nebu¬ 
lae: An Eye to the Euture, ed. Manchado A., Stanghellini 
L., Sch onberner D., (Gambridge:Gambridge Univ. Press), 
156 

Soker N., Livio M., 1994, ApJ, 421, 219 
Steffen W., Lopez J. A., Lim A., 2001, ApJ, 556, 823 
Tocknell J., De Marco O., Wardle M., 2014, MNRAS, 439, 
2014 

Tylenda R., Acker A., Stenholm B., Koeppen J., 1992, 
AVAS, 95, 337 

Tylenda R., Acker A., Stenholm B., 1993, AVAS, 102, 595 
Wang W., Liu X.-W., Zhang Y., Barlow M. J., 2004, AVA, 
427, 873 

Wesson R., Liu X.-W., Barlow M. J., 2005, MNRAS, 362, 
424 


© 0000 RAS, MNRAS 000, [^T? 


Low-ionization structures in planetary nebulae 


19 


Table 2. IC 4846 


Line ID 

Inner NEB 

NE jet 

SW jet 

Entire NEB 


3.5"(5)“ 

2.l"(3) 

2.l"(3) 

10"(15) 

[O II] 3727 

20.45 

100.2 

116.5 

26.78 

H12 3750 

2.421 

— 

— 

2.682 

Hll 3770 

2.994 

— 

— 

3.364 

HIO 3798 

4.415 

— 

— 

4.723 

He I 3819 

0.933 

— 

— 

1.216 

H9 3835 

5.792 

— 

— 

6.461 

[Ne III] 3869 

66.91 

84.11 

102.9 

76.43 

[Ne III] + Hi 3968 

36.53 

40.65 

56.82 

40.52 

He I 4026 

2.086 

2.986 

3.812 

2.403 

[S II] 4070 

1.553 

3.014 

2.191 

1.655 

H(5 4101 

24.41 

29.33 

32.56 

27.51 

H7 4340 

46.83 

50.58 

52.54 

52.34 

[O III] 4363 

7.274 

7.952 

7.768 

7.866 

He I 4471 

4.152 

6.811 

6.203 

4.885 

He II 4686 

0.461 

— 

— 

0.526 

He I 4712“ 

0.606 

— 

0.905 

0.713 

[Ar IV] 4712 

1.035 

— 

3.026 

1.111 

[Ar IV] 4740 

1.349 

— 

2.587 

1.453 

H/3 4861 

100.0 

100.0 

100.0 

100.0 

He I 4924 

1.198 

— 

— 

1.267 

[O III] 4958 

378.2 

349.5 

337.5 

421.1 

[O III] 5007 

1150 

1116 

982.9 

1257 

[N i] 5200 

0.173 

— 

— 

0.201 

[Fe III] 5270 

0.129 

— 

— 

0.174 

[Cl III] 5517 

0.333 

0.852 

1.328 

0.377 

[Cl III] 5537 

0.466 

0.771 

1.227 

0.546 

[N II] 5755 

0.425 

1.282 

1.089 

0.459 

He I 5876 

14.19 

15.59 

14.14 

14.18 

[K iv] 6101 

0.093 

— 

— 

0.108 

[O i] 6300 

1.676 

6.549 

5.691 

1.993 

[S III] 6312 

1.141 

1.512 

1.646 

1.228 

[O i] 6363 

0.582 

2.178 

2.238 

0.664 

[N II] 6548 

5.159 

20.79 

17.59 

6.487 

Ha 6564 

285.0 

285.0 

285.0 

285.0 

[N II] 6584 

16.54 

51.45 

45.97 

19.13 

He I 6678 

3.142 

3.546 

3.769 

3.089 

[S II] 6716 

1.128 

5.433 

4.996 

1.378 

[S II] 6731 

2.069 

7.506 

7.052 

2.579 

log[F(H/3) 

-11.62 

-13.49 

-13.45 

-11.54 

c(H/3) 

0.48±0.03 

0.50±0.05 

0.53±0.05 

0.52±0.03 

Line Fluxes 


Percentage errors in line fluxes 


(0.001-0.01)Fh^ 

13 

43 

43 

12 

(0.01-0.05)Fh^ 

10 

27 

26 

10 

(0.05-0.15)Fjj^ 

9 

12 

12 

8 

(0.15-0.30)Fh/3 

8 

10 

10 

7 

(0.30-2.0)Fh^ 

7 

9 

9 

7 

(2.0-5.0)FHy3 

7 

8 

7 

7 

(5.0-10.0)Fh^ 

6 

7 

7 

6 

> 10 

6 

7 

7 

6 


Electron 

Densities (cm~ 

and Temperatu 

res (K) 

Ne [S II] 

7200±1050 

2150±400 

2350±450 

8050±1050 

A/'e [Ar iv] 

8550±1400 

— 

2200±800 

9100±1500 

Ae[Cl III] 

7400±1350 

1650±950 

1850±700 

8050±1500 

Te [O III] 

9900±1400 

10350±1800 

11100±1950 

9850±1400 

Te [N II] 

11950±2250 

12800±4150 

12450±4000 

11400±2150 

Te[S II] 

12700±2000 

13700±4300 

13250±4150 

9500±1550 


^ The size of the component in arcsec. The number of pixels used to extract the spectra is indicated in the parentheses for each component. 


Table 3. Ionic and total abundances of IC 4846 


Ionic abundances 

Te & at « 

Inner NEB 

NE jet 

SW jet 

Entire NEB 

He+/H + 

1 

0.077±0.009 

0.091±0.014 

0.088±0.013 

0.078±0.009 

He++ /H + 

2 

0.0004±0.0001 

- 

- 

0.0004±0.0001 

ICF(He) 


1.00 

1.00 

1.00 

1.00 

He/H 


0.077±0.009 

0.091±0.014 

0.088±0.013 

0.078±0.009 

10^ X 0+/H + 

1 

0.78±0.08 

2.07±0.28 

2.55±0.33 

1.31±0.13 

10 ^ X o2+/h+ 

2 

40.7±3.8 

34.4±3.9 

23.9±2.8 

45.4±4.2 

ICF(O) 


l.OlitO.Ol 

1.00 

1.00 

l.OlitO.Ol 

10^ X O/H 


4.16±0.40 

3.65±0.42 

2.65±0.31 

4.69±0.43 

10^ X N+ /H + 

1 

0.22±0.03 

0.61±0.08 

0.56±0.07 

0.29±0.03 

ICF(N) 


53.3±7.6 

17.6±3.1 

10.4±1.8 

35.8±4.9 

10^ X N/H 


1.17±0.24 

1.06±0.24 

0.58±0.13 

1.04±0.18 

10® X S+/H + 

1 

0 .12±0.02 

0.25±0.07 

0.24±0.07 

0.16±0.02 

10 ® X s2+/h+ 

2 

2.75±0.26 

3.18±0.97 

2.53±0.71 

2.98±0.36 

ICF(S) 


2.63±0.12 

1.84±0.10 

1.56±0.08 

2.31±0.10 

10® X S/H 


7.53±0.82 

6.31±1.95 

4.33±1.24 

7.24±0.93 

10® X Cl2+ /H + 

2 

6 .11±0.86 

10.1±4.3 

10.1±3.8 

7.25±1.07 

ICF(Cl) 


2.74±0.40 

1.98±0.86 

1.71±0.69 

2.43±0.43 

lo"^ X Cl/H 


1.68±0.34 

2 .00±1.22 

1.73±0.95 

1.76±0.41 

lo"^ X Ar®+ /H + 

2 

2.35±0.33 

- 

3.75±1.16 

2.55±0.35 

ICF(Ar) 


1 .02±0.20 

1.06±0.26 

1.11±0.27 

1.03±0.20 

lo"^ X Ar/H 


2.40±0.59 

- 

4.15±1.64 

2.63±0.57 

10® X Ne2+/H + 

2 

7.03±0.69 

7.51±0.96 

6.89±0.88 

8.14±0.81 

ICF(Ne) 


1.02±0.14 

1.06±0.17 

1.11±0.18 

1.03±0.14 

10® X Ne/H 


7.19±1.21 

7.96±1.65 

7.62±1.59 

8.40±1.38 


“ 1: Te[N II] and ATg [S II], 2: Tg [O III] and Afg [S II] 
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Table 4. Wray 17-1 (PA=155°) 


Line ID 

NW filament 
11.6"(29) 

NW knot 
4.4"(11) 

NW inner NEB 
11.6"(29) 

SE filament 
10.8"(27) 

SE knot 
3.6" (9) 

SE inner NEB 
11.6"(29) 

[O II] 3727 

611.3 

364.2 

_ 

636.1 

309.3 

_ 

[Ne in] 3869 

137.7 

107.7 

— 

160.8 

112.5 

— 

[Ne III] + Hi 3968 

54.22 

— 

— 

65.39 

43.71 

— 

H(5 4101 

25.57 

— 

— 

25.57 

— 

— 

Hq 4340 

46.09 

48.56 

45.31 

49.85 

— 

44.21 

[O III] 4363 

14.18 

10.48 

6.044 

16.37 

— 

8.292 

He I 4471 

2.666 

— 

— 

1.587 

— 

— 

He II 4686 

53.93 

66.46 

113.6 

60.58 

63.21 

112.8 

[Ar iv]+He i 4712 

2.795 

— 

4.921 

2.456 

— 

4.056 

[Ar IV] 4712 

2.516 

— 

— 

2.289 

— 

— 

[Ar IV] 4740 

1.817 

— 

3.505 

1.659 

— 

2.845 

H/3 4861 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

He I 4924 

1.009 

— 

— 

1.016 

— 

— 

[O III] 4958 

405.3 

298.4 

118.1 

466.5 

393.4 

154.9 

[O III] 5007 

1174 

927.9 

354.1 

1369 

1191 

470.2 

[N i] 5200 

6.388 

— 

— 

5.296 

— 

— 

He II 5411 

4.601 

— 

— 

3.846 

— 

— 

[N II] 5755 

7.042 

— 

— 

6.536 

— 

— 

He I 5876 

9.481 

— 

— 

9.357 

— 

— 

[O i] 6300 

40.43 

63.24 

— 

32.23 

43.38 

— 

[O i] 6363 

14.66 

17.32 

— 

10.34 

18.43 

— 

[N II] 6548 

112.8 

78.61 

— 

92.77 

— 

— 

Ho; 6564 

285.0 

285.0 

285.0 

285.0 

285.0 

285.0 

[N II] 6584 

358.1 

241.8 

- 

300.9 

119.8 

- 

He I 6678 

2.545 

— 

— 

3.031 

— 

— 

[S II] 6716 

17.82 

6.849 

— 

15.01 

— 

— 

[S II] 6731 

15.71 

5.911 

— 

12.89 

— 

— 

He I 7065 

2.884 

— 

— 

2.493 

— 

— 

[Ar III] 7135 

24.26 

19.71 

— 

24.64 

19.88 

— 

[O II] 7330 

24.11 

— 

— 

— 

— 

— 

[Ar III] 7752 

3.396 

- 

- 

4.081 

- 

- 

log[T(H/3)] 

-13.52 

-15.05 

-14.51 

-13.52 

-15.22 

-14.26 

c(H/3) 

0.09±0.03 

0.10±0.08 

0.06±0.05 

0.11±0.04 

0.12±0.09 

0.07±0.06 

Line Fluxes 
(0.01-0.05 )Fh/3 

26 

_ 

Percentage errors in line fluxes 

45 25 

_ 

44 

(0.05-0.15)FHy3 

18 

43 

36 

17 

45 

33 

(0.15-0.30)Fjjy3 

12 

31 

26 

12 

35 

22 

(0.30-2.0)FHy3 

8 

18 

15 

9 

20 

15 

(2.0-5.0)Fh^ 

7 

16 

12 

8 

17 

11 

(5.0-10.0)Fh^ 

6 

14 

- 

6 

15 

- 

> 10 l^H/3 

6 

- 

- 

6 

12 

- 

Ne [S II] 

350±70 

Electron 

300±180 

Densities (cm~^ 

) and Temperatures (K) 
306±70 - 


Ne [Ar iv] 

300±110 

— 

160±110°' 

325±115 

— 

100±60°' 

Tg [O III] 

12400±2560 

12100±5490 

14250±5460 

12370±2460 

— 

14471±5250 

Tg [N II] 

11750±2350 

- 

- 

12360±2110 

- 

- 


^ Afg [Ar IV] is considered as a lower limit since the contribution of He I 4712 has not been subtracted. 


Table 5. Ionic and total abundances of Wray 17-1 (PA=155°) 


Ionic abundances 

Tg & A ® 

NW filament 

NW knot 

NW inner NEB 

SE filament 

SE knot® 

SE inner NEB 

He+/H+ 

1 

0.068±0.014 

_ 

_ 

0.071±0.014 

_ 

_ 

He++/H + 

2 

0.046±0.005 

0.056±0.013 

0.098±0.02 

0.047±0.006 

0.053±0.015 

0.097±0.021 

ICF(He) 


1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

He/H 


0.114±0.019 

0.056±0.013 

0.098±0.02 

0.118±0.020 

0.053±0.015 

0.097 

10^ X 0+/H + 

1 

12.5±1.3 

6.61±1.61 

- 

11.6±1.3 

5.98±1.57 

- 

10 ^ X o 2 +/h+ 

2 

20.7±2.1 

17.7±4.1 

4.21±0.82 

24.2±2.6 

23.4±5.6 

5.35±0.99 

ICF(O) 


1.41±0.09 

- 

- 

1.40±0.09 

- 

1.00 

10^ X O/H 


4.68±0.55 

>2.43±0.81 

>0.42±0.82 

5.02±0.63 

>2.93±1.04 

0.54±0.99 

10^ X N+/H + 

1 

4.53±0.48 

2.88±0.71 

- 

3.49±0.42 

1.48±0.42 

- 

ICF(N) 


3.76±0.58 

- 

- 

4.33±0.72 

- 

- 

10^ X N/H 


1.70±0.32 

>0.29±0.07 

- 

1.51±0.31 

>0.15±0.04 

- 

10® X S+/H + 

1 

0.56±0.09 

0 . 20 ± 0.10 

- 

0.41±0.07 

- 

- 

10 ® X s 2 +/h+ 

2 

3.01±0.48^ 

- 

- 

2.13±0.38^ 

- 

- 

ICF(S) 


1.18±0.05 

- 

- 

1.22±0.05 

- 

- 

10® X S/H 


4.23±0.68 

> 0 . 20 ± 0.10 

- 

3.21±0.57 

- 

- 

lo"^ X Ar2+/H + 

1 

12.8±1.9 

11.8±4.2 

- 

12 . 1 ± 1.8 

12.3±4.9 

- 

lo”^ X Ar®+/H + 

2 

2.12±0.58 

- 

2.94±1.39 

1.95±0.52 

- 

2.33±1.08 

ICF(Ar) 


1.36±0.11 

- 

- 

1.30±0.09 

- 

- 

10® X Ar/H 


2.04±0.38 

>1.18±0.42 

>0.29±0.14 

1.83±0.33 

> 1.23±0.50 

>0.23±0.11 

10® X Ne2+ /H + 

2 

6.17±0.69 

5.28±1.34 

- 

7.28±0.93 

5.86±1.18 

- 

ICF(Ne) 


2.26±0.35 

- 

- 

2.07±0.34 

- 

- 

10^ X Ne/H 


1.39±0.27 

>0.53±0.13 

- 

1.51±0.32 

>0.59±0.12 

- 


“ 1: Te[N II] and ATg [S II], 2: Tg [O III] and Alg [N II] 

^ Tg and Nq are drawn from the NW knot. 

^ It is calculated by using the equation =s+* ^4 0 yy_j_(-Q + 2 ^q^O. 433^ from Kingsburgh & Barlow (1994). 
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Table 6. Wray 17-1 (PA=75°) 


Line ID 

E knot 

2.8" (7) 

E arc 

14.8" (37) 

E inner NEB 
11.6" (29) 

W knot 

3.6" (9) 

W arc 
6"(15) 

W inner NEB 
11.6"(29) 

[O II] 3727 

237.4 

57.55 

_ 

253.3 

38.46 

_ 

[Ne III] 3869 

66.42 

52.69 

— 

113.9 

58.18 

— 

[Ne III] + Hi 3968 

28.35 

31.61 

— 

47.23 

21.86 

— 

H(5 4101 

24.98 

22.07 

— 

26.67 

23.69 

— 

Hy 4340 

46.89 

46.45 

44.89 

45.39 

47.89 

48.79 

[O III] 4363 

9.577 

9.244 

— 

9.039 

12.83 

4.863 

He II 4686 

92.72 

106.1 

103.8 

96.99 

100.1 

113.6 

[Ar lv]+He I 4712 

6.617 

6.087 

9.894 

7.862 

8.816 

10.53 

[Ar IV] 4712 

6.351 

— 

— 

7.691 

8.671 

— 

[Ar IV] 4740 

4.597 

4.366 

7.071 

5.668 

6.238 

7.414 

H/3 4861 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

He I 4924 

— 

— 

— 

5.043 

— 

— 

[O III] 4958 

180.9 

281.5 

85.08 

250.2 

308.9 

86.51 

[O III] 5007 

533.1 

843.8 

259.2 

762.5 

937.1 

263.1 

[N i] 5200 

4.978 

— 

— 

2.563 

— 

— 

He II 5411 

9.557 

8.051 

— 

7.658 

11.34 

3.934 

[N II] 5755 

3.842 

— 

— 

3.621 

— 

— 

He I 5876 

6.556 

— 

— 

3.958 

3.409 

— 

[O i] 6300 

37.13 

— 

— 

24.06 

— 

— 

[S III] 6312 

- 

- 

- 

2.541 

- 

- 

[O i] 6363 

11.05 

— 

— 

6.103 

— 

— 

[N II] 6548 

60.38 

8.051 

— 

49.81 

7.555 

— 

Ha 6564 

285.0 

285.0 

285.0 

285.0 

285.0 

285.0 

[N II] 6584 

183.2 

21.82 

- 

159.5 

26.55 

4.999 

He I 6678 

— 

— 

— 

— 

2.094 

— 

[S II] 6716 

4.986 

— 

— 

11.91 

5.129 

— 

[S II] 6731 

4.435 

— 

— 

12.36 

4.083 

— 

[Ar v] 7005 

5.591 

— 

7.921 

— 

— 

7.619 

[Ar III] 7135 

12.29 

12.58 

— 

15.42 

15.99 

3.479 

[O II] 7330 

- 

- 

- 

5.383 

- 

- 

log[F(H/3)] 

-14.85 

-14.36 

-14.49 

-14.51 

-14.41 

-14.29 

c(H/3) 

0.13±0.10 

0.09±0.07 

0.09±0.07 

0.11±0.06 

0.06±0.05 

0.04±0.04 

Line Fluxes 
(0.01-0.05 )Fh/3 

43 

38 

Percentage errors in line fluxes 
- 37 

36 

36 

(0.05-0.15)Fh^ 

36 

31 

34 

30 

28 

30 

(0.15-0.30)FHy3 

26 

26 

26 

20 

19 

23 

(0.30-2.0)Fh^ 

17 

14 

14 

12 

11 

14 

(2.0-5.0)Fh/3 

14 

10 

11 

9 

9 

11 

(5.0-10.0)Fh^ 

10 

9 

- 

8 

8 

- 

Ae [S II] 

375±230 

Electron 

Densities (cm 

®) and Temperatures (K) 
760±320 175±80 


Ae [Ar iv] 

325±200 

220 ± 110 ‘^ 

170±85°' 

515±220 

260±100 

100±45“ 

Te [O III] 

14590±5550 

12180±4025 

— 

12300±3920 

13030±3920 

14700±5600 

Te [N II] 

12130±5600 



12420±4880 




ATe [Ar IV] is considered as a lower limit since the contribution of He I 4712 has not been subtracted. 


Table 7. Ionic and total abundances of Wray 17-1 (PA=75°) 


Ionic abundances Tq 8z E knot E arc E inner NEB^ W knot W arc W inner NEB 


He+/H+ 1 

He++/H+ 2 


ICF(He) 

He/H 

10^ X 0+/H+ 1 

10^ X 0^+ /H+ 2 

ICF(O) 

10"^ X O/H 

10^ X N+ /H+ 1 

ICF(N) 

10"^ X N/H 

10® X S+/H+ 1 

10® X s2+/H+ 2 

ICF(S) 

10® X S/H 

lo"^ X Ar2+/H+ 2 

lo"^ X Ar®+/H+ 2 

lo"^ X Ar4+/H+ 2 

ICF(Ar) 

10® X Ar/H 

10® X Ne2+/H+ 2 

ICF(Ne) 

10^ X Ne/H 


0.047±0.020 

- 

0.083±0.020 

0.089±0.018 

1.000 

1.000 

0.130±0.040 

0.089±0.018 

4.34±0.97 

0.99±0.20 

6.01±1.19 

16.7±2.8 

1.97±0.41 

- 

2.04±0.60 

>1.77±0.46 

2.19±0.53 

0.26±0.08 

4.70±1.74 

— 

1.03±0.46 

>0.03±0.01 

0.16±0.08 

- 

0.84±0.41‘= 

- 

1.25±0.12 

- 

1.26±0.64 

- 

5.15±2.33 

8.33±2.83 

3.72±1.69 

5.37±2.18 

4.42±2.01 

- 

1.27±0.17 

- 

1.69±0.80 

>1.37±0.73 

1.79±0.43 

2.66±0.52 

3.39±1.21 

- 

0.61±0.26 

>0.27±0.05 


- 

0.029±0.011 

0.090±0.017 

0.084±0.014 

1.000 

1.000 

0.090±0.017 

0.113±0.025 

- 

4.52±0.68 

2.81±0.55 

13.6±1.9 

- 

2.47±0.51 

>0.28±0.06 

4.47±1.13 

- 

1.73±0.31 

- 

9.89±2.90 

- 

1.71±0.59 

- 

0.38 ±0.12 

- 

2.71±1.06 

- 

1.54±0.14 

- 

4.76±1.87 

- 

8.89±2.87 

5.66±2.13 

6.62±2.38 

6.15±2.31 

- 

- 

1.11±0.05 

>1.28±0.68 

>1.73±0.59 

- 

5.26±0.91 

- 

3.30±0.96 

- 

1.74±0.59 


0.032±0.012 

- 

0.085±0.013 

0.097±0.017 

1.000 

1.000 

0.117±0.025 

0.097±0.017 

0.45±0.08 

- 

14.2±1.9 

2.86±0.51 

2.37±0.47 

- 

3.47±0.83 

>0.29±0.05 

0.23±0.07 

0.04±0.02 

76.9±22.9 

- 

1.77±0.75 

>0.004±0.002 

0.12±0.05 

- 

0.64±0.27'= 

- 

2.96±0.29 

- 

2.26±0.97 

- 

8.27±2.11 

1.44±0.56 

6.51±2.02 

5.95±1.97 

- 

5.92±1.96 

1.01±0.07 

- 

>1.50±0.42 

>1.33±0.81 

2.23±0.35 

- 

2.45±0.67 

- 

0.55±0.17 

- 


“ l:re[N II] and Ag [S II], 2: Tg [O III] and Ag [N II] 

® Tq is drawn from the W nebula 

^ It is calculated by using the equation =s +*^4 0 yy_j_(-Q + 2 ^q^O. 433^ from Kingsburgh & Barlow (1994). 
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Table 8. K 1-2 


Line ID 

NW filament 

NW knot 

NW filament-|-Knot 

SE knot 

NW NEB 


7.6"(19) 

3.6" (9) 

11.6"(29) 

3.6" (9) 

35.6" (89) 

[O II] 3727 

272.1 

193.1 

245.2 

265.6 

205.3 

[Ne III] 3869 

117.1 

41.47 

94.17 

— 

83.57 

[Ne III] + Hi 3968 

57.85 

20.94 

48.83 

— 

38.51 

H5 4101 

23.63 

26.61 

24.31 

— 

23.93 

H 7 4340 

45.86 

50.61 

45.82 

49.21 

48.55 

[O III] 4363 

16.87 

7.828 

12.95 

7.118 

13.94 

He II 4686 

75.32 

94.69 

79.09 

107.3 

89.62 

He I 4712" 

0.358 

— 

0.316 

0.509 

0.292 

[Ar IV] 4712 

6.544 

8.405 

8.430 

9.265 

7.599 

[Ar IV] 4740 

4.916 

— 

4.955 

6.801 

3.673 

H/3 4861 

100.0 

100.0 

100.0 

100.0 

100.0 

[O III] 4958 

296.9 

155.3 

266.1 

146.6 

262.4 

[O III] 5007 

905.1 

475.1 

819.3 

417.4 

789.3 

[N i] 5200 

5.733 

— 

4.049 

— 

4.948 

He II 5411 

5.922 

— 

4.362 

— 

6.499 

[N II] 5755 

2.919 

— 

2.628 

— 

2.759 

He I 5876 

8.443 

— 

6.523 

— 

6.872 

[O i] 6300 

22.84 

22.38 

23.14 

20.28 

15.23 

[S III] 6312 

8.126 

7.638 

10.24 

— 

10.49 

[O i] 6363 

7.581 

10.88 

8.006 

— 

5.313 

[N ii] 6548 

87.19 

57.08 

85.35 

49.29 

62.24 

Ha 6564 

285.0 

285.0 

285.0 

285.0 

285.0 

[N ii] 6584 

270.1 

184.1 

263.3 

145.7 

208.9 

He i 6678 

2.938 

— 

3.476 

— 

— 

[S ii] 6716 

36.68 

19.91 

30.69 

16.26 

23.03 

[S ii] 6731 

35.47 

21.51 

30.99 

14.68 

23.81 

[Ar v] 7005 

2.687 

— 

2.239 

— 

4.491 

[Ar III] 7135 

16.22 

11.51 

15.59 

12.27 

15.52 

log[A'(H/3)] 

-14.21 

-14.87 

-14.11 

-15.01 

-13.94 

c(H/3) 

0.26±0.06 

0.25±0.07 

0.22±0.06 

0.25±0.08 

0.25±0.06 

Line Fluxes 


Percentage errors in line fluxes 


(0.01-0.05)Fh^ 

42 

51 

40 

55 

35 

(0.05-0.15)Ffj^ 

32 

35 

29 

38 

26 

(0.15-0.30)Fh/3 

16 

22 

15 

25 

14 

(0.30-2.0)Fh^ 

11 

14 

11 

16 

10 

(2.0-5.0)Fh;3 

9 

10 

9 

10 

8 

(5.0-10.0)Fh^ 

8 

- 

8 

- 

7 



Electron Densities (cm and Temperatures (K) 


Ne [S II] 

520±80 

888±280 

610±100 

422±150 

700±140 

Ne [Ar iv] 

755±340 

— 

— 

492±260 

— 

Te[0 III] 

14820±2880 

14010±5500 

13800±4300 

14250±5690 

14470±2500 

Te [N II] 

9120±3950 


8850±3700 


9820±3600 


Table 9. Ionic and total abundances of K 1-2 


Ionic abundances 

Te & AT “ 

NW filament 

NW knot 

NW filament-|-Knot 

SE knot 

NW NEB 

He+/H+ 

1 

0.068±0.023 

_ 

0.054±0.017 

_ 

0.049±0.015 

He + +/H + 

2 

0.068±0.011 

0.082±0.014 

0.067±0.010 

0.092±0.020 

0.077±0.011 

ICF(He) 


1.000 

1.00 

1.000 

1.000 

1.000 

He/H 


0.136±0.034 

0.082±0.014 

0.121±0.027 

0.092±0.020 

0.126±0.026 

10^ X 0+/H + 

1 

16.8±2.4 

2.35±0.46 

17.9±2.5 

2.81 ±0.43 

9.25±1.19 

10 ® X o 2 +/h+ 

2 

9.85±1.34 

5.83±1.02 

10.5±1.4 

5.07 ±0.96 

9.03±1.09 

ICF(O) 


1.59±0.20 

- 

1.71±0.22 

- 

1.88±0.26 

10^ X O/H 


4.23±0.79 

> 0.82±0.16 

4.86±0.92 

>0.79±0.19 

3.43±0.64 

10® X N+ /H + 

1 

6.64±1.08 

1.63±0.28 

7.07±1.01 

1.25±0.28 

4.11±0.55 

ICF(N) 


2.52±0.59 

- 

2.71±0.64 

- 

3.71±0.84 

10^ X N/H 


1.67±0.48 

>0.16±0.03 

1.92±0.53 

>0.13±0.03 

1.52±0.40 

10® X S+/H + 

1 

2.28±0.41 

0.53±0.15 

2.16±0.36 

- 

1.26±0.22 

10 ® X s 2 +/h+ 

2 

4.59±1.55 

5.18±1.95 

6.48±1.94 

- 

4.21±1.26 

ICF(S) 


1.09±0.05 

- 

1.10±0.05 

- 

1.18±0.07 

10® X S/H 


7.46±2.15 

>5.71±1.26 

9.52±2.57 

- 

6.45±1.78 

lo"^ X Ar2+/H + 

2 

6.59±1.49 

5.21±1.96 

7.24±2.11 

5.38±2.21 

6.61±1.98 

lo"^ X Ar®+ /H + 

2 

3.78±1.56 

5.54±2.61 

5.08±1.97 

5.76±2.37 

3.71±1.23 

lo"^ X Ar4+/H + 

2 

2.05±0.86 

- 

2.01±0.83 

- 

3.61±1.37 

ICF(Ar) 


1.66±0.36 

- 

1.59±0.29 

- 

1.37±0.15 

10® X Ar/H 


2.06±0.79 

>1.08±0.65 

2.27±0.88 

>1.11±0.65 

1.91±0.66 

10® X Ne2+/H + 

2 

3.01±0.47 

1.26±0.27 

3.01±0.47 

- 

2.31±0.41 

ICF(Ne) 


4.30±0.99 

- 

4.64±1.08 

- 

3.80±0.84 

10^ X Ne/H 


1.29±0.48 

>0.13±0.03 

1.40±0.39 

- 

0.88±0.25 


“ 1: Te[N II] and ATg [S II], 2: Tg [O III] and Afg [S II] 
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Table 10. NGC 6891 (PA=135°) 


Line ID 

Inner NEB 
4 . 9 " (7) 

NW knot 
2.l"(3) 

SE knot 
2.l"(3) 

NW outer NEB 
4 . 9 " (7) 

SE outer NEB 
4.9"(7) 

Entire NEB 
20 . 3 " (29) 

[O II] 3727 

12.89 

34.87 

24.18 

22.59 

18.15 

16.55 

H12 3750 

— 

— 

— 

3.345 

2.239 

2.266 

Hll 3770 

— 

— 

— 

4.037 

3.069 

2.978 

HIO 3798 

— 

— 

— 

5.341 

4.001 

3.381 

He I 3819 

— 

— 

— 

1.191 

1.161 

1.489 

H9 3835 

— 

— 

— 

7.707 

5.777 

4.999 

[Ne III] 3869 

55.17 

62.19 

45.15 

56.17 

43.27 

52.34 

[Ne III] + Hi 3968 

26.02 

32.13 

26.71 

31.89 

26.73 

28.41 

He I 4026 

1.886 

— 

— 

2.129 

2.013 

1.876 

[S ii] 4070 

1.184 

— 

— 

— 

— 

0.721 

H5 4101 

22.21 

24.81 

26.64 

24.92 

24.69 

23.55 

Hq 4340 

43.62 

46.12 

49.89 

50.29 

48.47 

44.14 

[O III] 4363 

5.311 

4.811 

4.669 

4.894 

4.855 

4.817 

He I 4387 

0.723 

— 

— 

1.088 

0.639 

0.609 

He I 4471 

4.755 

6.338 

5.978 

5.856 

4.341 

4.973 

He II 4686 

1.817 

— 

— 

— 

— 

0.631 

He I 4712*^ 

0.499 

0.665 

0.628 

0.615 

0.465 

0.522 

[Ar iv] 4712 

0.625 

— 

— 

0.455 

0.508 

0.597 

[Ar iv] 4740 

0.756 

— 

— 

0.439 

0.419 

0.518 

H/3 4861 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

He I 4924 

1.055 

— 

— 

1.441 

1.302 

1.322 

[O III] 4958 

296.2 

280.9 

296.6 

266.5 

290.6 

289.1 

[O III] 5007 

851.3 

833.9 

885.3 

813.4 

893.6 

858.6 

[Cl III] 5517 

0.232 

— 

— 

0.437 

0.357 

0.488 

[Cl III] 5537 

0.223 

— 

— 

0.431 

0.363 

0.501 

[N II] 5755 

— 

— 

— 

0.107 

0.091 

0.167 

He I 5876 

16.86 

14.67 

17.74 

15.09 

17.44 

16.38 

[S III] 6312 

0.292 

— 

— 

0.623 

0.509 

0.361 

[N II] 6548 

2.124 

6.153 

5.272 

2.029 

2.807 

2.752 

Ha 6564 

285.0 

285.0 

285.0 

285.0 

285.0 

285.0 

[N II] 6584 

5.771 

14.88 

16.92 

7.335 

7.684 

6.971 

He I 6678 

4.479 

4.883 

5.386 

4.151 

3.985 

4.333 

[S II] 6716 

0.261 

0.799 

1.473 

0.229 

0.355 

0.271 

[S II] 6731 

0.312 

0.929 

1.559 

0.314 

0.496 

0.332 

log[F(H/3)] 

-11.68 

-13.58 

-13.52 

-12.14 

-12.05 

-11.29 

c(H/3) 

0.03±0.02 

0.03±0.03 

0.03±0.03 

0.06±0.04 

0.05±0.03 

0.05±0.03 

Line Fluxes 

(0.001-0.01)Ffj^ 

32 

Percentage errors in line fluxes 

48 49 36 

36 

30 

(0.01-0.05)Ffj/3 

17 

36 

35 

22 

23 

15 

(0.05-0.15)Fjj^ 

13 

32 

30 

16 

17 

11 

(0.15-0.30)Ffj^ 

9 

22 

21 

13 

12 

9 

(0.30-2.0)FHy3 

8 

13 

14 

9 

8 

7 

(2.0-5.0)Ffjy3 

6 

11 

10 

8 

8 

6 

Ne [S II] 

Electron 

1190±530 

Densities (cm 
1150±670 

and Temperatures (K) 
750±440 1910±980 

1980±1020 

1420±600 

Ne [Ar iv] 

7470±3400 

— 

— 

3650±1860 

1680±860 

2260±960 

Ag[Cl III] 

2180±980 

— 

— 

2430±1240 

2660±1360 

2820±1200 

Tg [O III] 

9930±1960 

9700±3690 

9460±3500 

9820±2420 

9530±2340 

9600±1740 

Tg [N II] 




10055±4130 

9200±3790 

12720±4230 


Table 11. Ionic and total abundances of NGC 6891 (PA=135°) 


Ionic abundances 

Te & 

Inner NEB 

NW knot 

SE knot 

NW outer NEB 

SE outer NEB 

Entire NEB 

He+/H+ 

1 

0.116±0.018 

0.104±0.040 

0.129±0.050 

0.103±0.016 

0.114±0.018 

0.107±0.015 

He++/H + 

2 

0.0015±0.0003 

- 

- 

- 

- 

0.0005±0.0002 

ICF(He) 


1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

He/H 


0.118±0.018 

0.104±0.040 

0.129±0.050 

0.103±0.016 

0.114±0.018 

0.108±0.015 

10^ X 0+/H + 

1 

0.59±0.09 

1.79±0.54 

1.32±0.41 

1.09±0.18 

1.34±0.21 

0.31±0.05 

10^ X 0^+ /H + 

2 

30.3±3.1 

31.9±7.8 

37.1±9.1 

29.3±3.5 

35.8±4.3 

34.1±3.7 

ICF(O) 


l.OlitO.OO 

1.00 

1.00 

1.00 

1.00 

l.OOitO.OO 

10^ X O/H 


3.12±0.32 

3.36±0.84 

3.84±0.95 

3.03±0.37 

3.71±0.45 

3.45±0.38 

10^ X N+/H + 

1 

0 . 12 ± 0.02 

0.35±0.14 

0.39±0.15 

0.15±0.03 

0.20±0.04 

0.08±0.02 

ICF(N) 


53.8±9.7 

18.8±7.4 

29.1±11.6 

27.8±5.7 

27.7±8.1 

111±22 

10^ X N/H 


0.63±0.16 

0.66±0.37 

1.13±0.63 

0.42±0.12 

0.55±0.16 

0.89±0.28 

lo"^ X S+/H + 

1 

0.16±0.05 

0.51±0.26 

0.91±0.46 

0.17±0.04 

0.33±0.08 

0.11±0.03 

10® X s2+/H + 

2 

0.71±0.24 

0.29±0.15® 

0.52±0.26® 

1.59±0.61 

1.49±0.56 

1.02±0.32 

ICF(S) 


2.62±0.16 

1.88±0.23 

2.16±0.28 

2.13±0.14 

2.12±0.14 

3.34±0.22 

10® X S/H 


1.90±0.65 

0.64±0.33 

1.31±0.68 

3.42±1.33 

3.23±1.22 

3.45±1.03 

10^ X CP+/H + 

2 

2.96±0.99 

- 

- 

5.95±2.22 

5.47±2.03 

7.25±2.27 

ICF(Cl) 


2.68±1.29 

2.21±1.62 

2.54±1.85 

2.15±1.17 

2.17±1.16 

3.38±0.15 

10^ X Cl/H 


7.93±4.65 

- 

- 

12.79±8.45 

11.87±7.71 

24.47±13.39 

lo"^ X Ar®+ /H + 

2 

1.29±0.43 

- 

- 

0.86±0.33 

0.97±0.35 

1.15±0.38 

ICF(Ar) 


1.03±0.19 

1.06±0.59 

1.04±0.58 

1.04±0.26 

1.04±0.26 

1 . 01 ± 0.20 

lo"^ X Ar/H 


1.33±0.51 

- 

- 

0.89±0.41 

1.01±0.44 

1.16±0.45 

10® X Ne2+ /H + 

2 

5.67±0.65 

7.09±2.21 

5.74±1.79 

6.07±0.77 

5.33±0.69 

6.24±0.57 

ICF(Ne) 


1.03±0.15 

1.06±0.37 

1.04±0.36 

1.04±0.18 

1.04±0.18 

1.01±0.16 

10® X Ne/H 


5.82±0.98 

7.49±3.51 

5.94±2.79 

6.30±1.34 

5.53±1.19 

6.32±1.13 


^ 1: Te[N II] and ATg [S II], 2: Tg [O III] and Afg [S II]. 

^ It is calculated by using the equation =g +*(-4 0 yy_j _^q + 2y q^O. 433^ from Kingsburgh & Barlow (1994). 
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Table 12. NGC 6891 (PA=45°) 


Line ID 

Inner NEB 

NE NEB 

SW NEB 

Entire NEB 


4.9"(7) 

6.3" (9) 

6.3" (9) 

20.3"(29) 

[O II] 3727 

11.31 

15.92 

17.08 

12.77 

H12 3750 

— 

2.752 

2.714 

2.326 

Hll 3770 

— 

3.204 

3.144 

2.633 

HIO 3798 

— 

4.656 

4.713 

3.099 

He I 3819 

— 

1.978 

1.573 

1.259 

H9 3835 

— 

6.902 

7.545 

4.285 

[Ne III] 3869 

54.81 

45.61 

48.82 

52.55 

[Ne III] + Hi 3968 

24.85 

27.59 

28.06 

27.19 

He I 4026 

1.698 

2.389 

2.206 

2.004 

[S II] 4070 

— 

— 

— 

0.474 

HS 4101 

22.84 

24.89 

24.07 

22.98 

H-v 4340 

42.87 

46.56 

43.87 

45.32 

[O III] 4363 

4.931 

4.152 

3.698 

4.682 

He I 4387 

0.591 

0.909 

0.731 

0.399 

He I 4471 

5.176 

5.317 

4.309 

5.045 

He II 4686 

1.994 

— 

— 

0.855 

He I 4712°= 

0.543 

0.558 

0.452 

0.529 

[Ar IV] 4712 

0.589 

0.749 

0.421 

0.705 

[Ar iv] 4740 

0.699 

— 

0.399 

0.645 

H/3 4861 

100.0 

100.0 

100.0 

100.0 

He I 4924 

1.263 

1.515 

1.061 

1.253 

[O III] 4958 

289.3 

279.5 

265.6 

302.3 

[O III] 5007 

889.9 

849.5 

780.7 

872.6 

[Cl III] 5517 

0.211 

0.407 

0.411 

0.529 

[Cl III] 5537 

0.198 

0.379 

0.397 

0.518 

[N II] 5755 

— 

— 

— 

0.158 

He I 5876 

16.94 

17.36 

15.05 

16.57 

[S III] 6312 

0.235 

0.338 

0.538 

0.311 

[N II] 6548 

2.393 

2.923 

2.512 

2.848 

Hq! 6564 

285.0 

285.0 

285.0 

285.0 

[N II] 6584 

5.916 

6.901 

6.953 

6.844 

He I 6678 

4.291 

4.194 

4.141 

4.348 

[S II] 6716 

0.225 

0.291 

0.209 

0.302 

[S II] 6731 

0.277 

0.374 

0.248 

0.349 

log[F(H/3)] 

-11.64 

-12.31 

-12.11 

-11.31 

c(H/3) 

0.05±0.03 

0.04±0.03 

0.05±0.03 

0.04±0.03 

Line Fluxes 


Percentage errors in line fluxes 


(0.001-0.01)Fjj^ 

33 

35 

34 

29 

(0.01-0.05)Eh/3 

17 

24 

25 

14 

(0.05-0.15)Fjjy3 

14 

17 

17 

12 

(0.15-0.30)Fh/3 

9 

14 

13 

9 

(0.30-2.0)Eh^ 

7 

8 

8 

7 

(2.0-5.0)Eh/3 

7 

8 

8 

7 


Electron 

Densities (cm 

®) and Temperatures (K) 

Ae [S II] 

1300±610 

1500±750 

1110±530 

1200±500 

Ae [Ar iv] 

7000±3270 

— 

3330±1600 

2900±1190 

Ae [Cl III] 

1950±910 

1920±960 

2420±1160 

2350±970 

Te [O III] 

9590±1920 

9250±2440 

9200±2510 

9500±1670 

Te [N II] 




12190±3870 


Table 13. Ionic and total abundances of NGC 6891 (PA=45°) 


Ionic abundances Tq & Inner NEB NE NEB SW NEB Entire NEB 


He+/H+ 1 

He + +/H+ 2 

ICF(He) 

He/H 

10^ X 0+ /H+ 1 

10^ X 0^+ /H+ 2 

ICF(O) 

10^ X O/H 

10^ X N+ /H+ 1 

ICF(N) 

10^ X N/H 

lo"^ X S+/H+ 1 

10® X 3^+ /H+ 2 

ICF(S) 

10® X S/H 

10® X CP+/H+ 2 

ICF(Cl) 

10® X Cl/H 

lo"^ X Ar®+ /H+ 2 

ICF(Ar) 

lo"^ X Ar/H 

10® X Ne2+/H+ 2 

ICF(Ne) 

10® X Ne/H 


0.113±0.017 

0.104±0.019 

0.0016±0.0004 

- 

1.000 

1.000 

0.115±0.017 

0.104±0.019 

0.63±0.09 

1.05±0.20 

34.9±3.5 

39.1±4.4 

l.OlitO.OO 

1.00 

3.58±0.36 

4.02±0.46 

0.14±0.03 

0.19±0.04 

56.9±9.9 

38.3±8.5 

0.80±0.22 

0.73±0.22 

0.16±0.05 

0.23±0.08 

0.67±0.24 

1.16±0.41 

2.68±0.15 

2.36±0.17 

1.84±0.67 

2.79±1.01 

2.99±1.03 

6.58±2.31 

2.75±1.398 

2.40±1.13 

8.21±5.046 

15.81±9.73 

1.35±0.46 

1.57±0.54 

1.03±0.21 

1.03±0.23 

1.39±0.55 

1.61±0.66 

6.56ib0.65 

6.56±0.74 

1.03±0.14 

1.03±0.17 

6.74±1.16 

6.74±1.33 


0.099±0.018 

0.109±0.015 

- 

0.0007±0.0002 

1.000 

1.000 

0.099±0.018 

0.110±0.015 

1.12±0.21 

0.26±0.04 

36.6±4.1 

36.5±3.6 

1.00 

l.OOibO.OO 

3.77±0.44 

3.70±0.37 

0.18±0.04 

0.09±0.02 

33.7±7.4 

142±26 

0.61±0.19 

1.28±0.37 

0.16±0.06 

0.21±0.06 

1.87±0.66 

0.92±0.29 

2.26±0.16 

3.63±0.22 

4.27±1.54 

3.41±1.09 

6.71±2.34 

7.89±2.35 

2.28±1.15 

3.71±1.67 

15.31±9.38 

29.26±15.77 

0.98±0.33 

1.45±0.44 

1.03±0.24 

1.01±0.20 

1.01±0.41 

1.47±0.53 

7.03±0.80 

6.56±0.65 

1.03±0.17 

1.01±0.14 

7.25±1.43 

6.64±1.13 


“ 1: Te[N II] and Ag [S II], 2: Tg [O III] and Ag [S II] 
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Table 14. NGC 6572 (PA=160°) 


Line ID 

Inner NEB 
7.7"(11) 

N knot 
2.l"(3) 

S knot 
2.l"(3) 

Entire NEB 
17.5"(25) 

[O II] 3727 

27.25 

127.8 

143.6 

34.75 

H12 3750 

3.901 

— 

— 

3.325 

Hll 3770 

4.193 

— 

— 

3.815 

HIO 3798 

5.295 

— 

— 

5.671 

He I 3819 

1.735 

— 

— 

1.771 

H9 3835 

7.183 

— 

— 

7.059 

[Ne III] 3869 

93.02 

111.3 

137.3 

96.82 

[Ne III] + Hi 3968 

45.14 

54.09 

72.41 

46.21 

He I 4026 

2.562 

— 

— 

2.739 

[S II] 4070 

2.537 

— 

— 

2.754 

H(5 4101 

23.95 

24.43 

28.53 

26.42 

Hq 4340 

46.22 

44.51 

45.26 

47.79 

[O III] 4363 

9.183 

7.143 

9.579 

9.651 

He I 4387 

0.805 

— 

— 

0.823 

He I 4471 

5.541 

6.701 

7.006 

5.827 

He II 4686 

1.014 

— 

— 

0.906 

He I 4712°= 

0.808 

0.978 

1.022 

0.851 

[Ar IV] 4712 

1.259 

— 

— 

1.172 

[Ar IV] 4740 

2.761 

— 

— 

2.607 

H/3 4861 

100.0 

100.0 

100.0 

100.0 

He I 4924 

1.475 

— 

— 

1.355 

[O III] 4958 

383.2 

372.9 

441.9 

418.8 

[O III] 5007 

1225 

1122 

1331 

1243 

[N i] 5200 

0.274 

5.733 

— 

0.348 

[Cl III] 5517 

0.239 

— 

— 

0.268 

[Cl III] 5537 

0.538 

— 

— 

0.602 

[N II] 5755 

2.136 

4.921 

3.665 

2.337 

He I 5876 

18.34 

17.98 

18.79 

18.27 

[O i] 6300 

5.225 

11.96 

17.97 

5.963 

[S III] 6312 

1.174 

— 

— 

1.201 

[O i] 6363 

1.684 

4.483 

6.707 

1.841 

[N II] 6548 

18.41 

67.85 

54.21 

19.81 

Ha 6564 

285.0 

285.0 

285.0 

285.0 

[N II] 6584 

55.56 

196.1 

161.6 

60.55 

He I 6678 

4.535 

4.045 

4.283 

4.524 

[S II] 6716 

0.763 

3.545 

2.867 

0.887 

[S II] 6731 

1.592 

6.832 

5.579 

1.853 

log[B(H/3)] 

-10.32 

-12.49 

-12.54 

-10.26 

c(H^) 

0.27±0.01 

0.25±0.05 

0.29±0.06 

0.27±0.02 

Line Fluxes 
( 0 . 001 - 0 . 01 )Fh /3 

13 

Percentage errors in line fluxes 
40 40 

13 

(0.01-0.05)Fh^ 

9 

33 

32 

9 

(0.05-0.15)Fh/3 

7 

23 

22 

7 

(0.15-0.30)Fh^ 

6 

18 

18 

6 

(0.30-2.0)Fh^ 

5 

12 

12 

6 

(2.0-5.0)Fh;3 

5 

9 

9 

5 

(5.0-10.0)Fh^ 

5 

- 

- 

5 

> 10 l^H/3 

5 

6 

6 

5 


Electron Densities (cm 

and Temperatures (K) 

Ne [S II] 

21300±3370 

9820±3900 

10420±4200 

21580±3400 

Ne [Ar iv] 

27000±3630 

— 

— 

29300±3940 

Ae[Cl III] 

35230±6480 

— 

— 

35200±6480 

Te [O III] 

10200±1230 

9900±2670 

10260±2770 

10200±1230 

Te [N II] 

12350±1640 

11350±3970 

10730±3750 

12220±1630 

Te[S II] 

9100±1300 



9820±1420 


Table 15. Ionic and total abundances of NGC 6572 (PA=160°) 


Ionic abundances 

Te & V “ 

Inner NEB 

N knot 

S knot 

Entire NEB 

He+/H+ 

1 

0.095±0.008 

0.099±0.021 

0.105±0.023 

0.095±0.008 

He + +/H + 

2 

0.0008±0.0001 

- 

- 

0.0008±0.0001 

ICF(He) 


1.00 

1.00 

1.0 

1.00 

He/H 


0.096±0.008 

0.099±0.021 

0.105±0.029 

0.096±0.008 

10^ X 0+/H + 

1 

1.77±0.13 

7.02±1.19 

10.21±1.73 

2.38±0.20 

10^ X 0^+ /H + 

2 

38.7±2.7 

36.1±5.4 

41.7±6.3 

40.9±3.2 

ICF(O) 


l.OlitO.ll 

1.00 

1.00 

l.OlitO.ll 

10^ X O/H 


4.07±0.54 

4.31±0.66 

5.20±0.80 

4.35±0.59 

10^ X N+/H + 

1 

0.76±0.06 

3.12±0.53 

2.93±0.49 

0.88±0.07 

ICF(N) 


23.0±3.5 

6.14±1.40 

5.09±1.16 

18.3±2.9 

10^ X N/H 


1.75±0.30 

1.915±0.55 

1.49±0.42 

1.608±2.88 

10® X S+/H + 

1 

0.17±0.02 

0.48±0.17 

0.46±0.16 

0.18±0.02 

10® X S2+ /H + 

2 

2.38±0.27 

2.69±0.95® 

2.57±0.90® 

2.45±0.28 

ICF(S) 


2 . 00 ± 0.10 

1.34±0.09 

1.28±0.08 

1.86±0.09 

10® X S/H 


5.10±0.63 

4.25±1.53 

3.87±1.37 

4.89±0.61 

10® X CP+/H + 

3 

6.74±0.94 

- 

- 

7.57±1.08 

ICF(Cl) 


2.14±0.36 

1.58±0.80 

1.51±0.75 

1.20±0.34 

10® X Cl/H 


14.45±0.32 

- 

- 

15.12±0.34 

lo"^ X Ar®+/H + 

2 

3.86±0.44 

- 

- 

3.59±0.42 

ICF(Ar) 


1.05±0.19 

1.20±0.39 

1.25±0.40 

1.06±0.20 

lo"^ X Ar/H 


4.06±0.87 

- 

- 

3.82±0.86 

10® X Ne2+/H + 

2 

8.58±0.61 

8.27±1.41 

12.31±2.11 

8.93±0.69 

ICF(Ne) 


1.05±0.16 

1.19±0.26 

1.25±0.27 

1.06±0.17 

10® X Ne/H 


9.02±1.50 

9.88±2.71 

15.31±4.21 

9.50±1.67 


“ 1: Te[N II] and ATg [S II], 2: Tg [O III] and A^e [Ar IV] and 3: Tg [O III] and Alg [Cl III]. 

^ It is calculated by using the equation S”t2 =g+* ^4 0 yy_j_ ^q + 2 y q^O .433 ^ from Kingsburgh & Barlow (1994). 
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Table 16. NGC 6572 (PA=15°) 


Line ID 

Inner NEB 
7.7"(11) 

N knot 
2.l"(3) 

S knot 
2.l"(3) 

Entire NEB 
17 . 5 " (25) 

[O ii] 3727 

27.91 

139.9 

101.3 

34.69 

H12 3750 

3.337 

— 

— 

2.773 

Hll 3770 

3.667 

— 

— 

3.531 

HIO 3798 

4.948 

— 

— 

4.986 

He I 3819 

1.243 

— 

— 

1.208 

H9 3835 

7.024 

— 

— 

6.933 

[Ne III] 3869 

92.18 

131.1 

113.6 

92.17 

[Ne III] 3968 

45.01 

66.21 

57.23 

44.48 

He I 4026 

2.729 

— 

— 

2.265 

[S II] 4070 

2.202 

— 

— 

2.332 

H5 4101 

25.75 

24.73 

26.06 

25.45 

Hq 4340 

45.62 

50.88 

44.76 

45.57 

[O III] 4363 

8.615 

11.71 

8.384 

8.675 

He I 4471 

5.399 

6.976 

9.133 

5.508 

He II 4686 

0.995 

— 

— 

0.853 

He I 4712°= 

0.788 

1.018 

1.331 

0.805 

[Ar IV] 4712 

1.169 

— 

— 

1.011 

[Ar IV] 4740 

2.558 

— 

— 

2.341 

H/3 4861 

100.0 

100.0 

100.0 

100.0 

He I 4924 

1.456 

— 

— 

1.646 

[O III] 4958 

420.1 

434.1 

387.8 

393.8 

[O III] 5007 

1229 

1360 

1104 

1207 

[N i] 5200 

0.373 

— 

— 

0.443 

[Cl III] 5517 

0.278 

— 

— 

0.276 

[Cl III] 5537 

0.559 

— 

— 

0.575 

[N II] 5755 

2.141 

3.171 

4.152 

2.316 

He I 5876 

17.84 

14.91 

14.95 

17.88 

[O i] 6300 

5.218 

15.08 

12.93 

6.112 

[S III] 6312 

1.038 

— 

— 

1.075 

[O i] 6363 

1.638 

4.837 

4.422 

1.913 

[N II] 6548 

16.78 

66.71 

58.81 

21.21 

Hq! 6564 

285.0 

285.0 

285.0 

285.0 

[N II] 6584 

52.88 

199.6 

177.1 

64.52 

He I 6678 

4.441 

3.887 

4.295 

4.424 

[S II] 6716 

0.731 

2.425 

2.394 

0.866 

[S II] 6731 

1.519 

4.682 

4.767 

1.803 

log[F(H^)] 

-10.31 

-12.24 

-12.35 

-10.23 

c(H/3) 

0.29±0.02 

0.29±0.06 

0.24±0.05 

0.28±0.01 

Line Fluxes 
(0.001-0.01 )Fh/3 

14 

Percentage errors in line fluxes 
38 38 

13 

(0.01-0.05)Fh^ 

9 

33 

33 

9 

(0.05-0.15)FHy3 

7 

22 

23 

7 

(0.15-0.30)Fh^ 

6 

17 

18 

6 

(0.30-2.0)Fh^ 

5 

13 

13 

5 

(2.0-5.0)Fh^ 

5 

8 

8 

5 

(5.0-10.0)Fh^ 

5 

- 

- 

5 

> 10 Fh;3 

5 

6 

6 

5 


Electron Densities (cm 

'^) and Temperatures (K) 

Ne [S II] 

20380±2600 

9600±4480 

12620±5900 

20380±2600 

Ne [Ar iv] 

26380±3360 

— 

— 

29650±3780 

Ne [Cl III] 

22750±4500 

— 

— 

25930±5130 

Tg [O III] 

10000±1200 

10890±2650 

10400±2620 

10050±1200 

Tg [N II] 

12780±1700 

9420±3275 

10570±3670 

11925±1600 

Te [S II] 

10360±1500 



8990±1300 


Table 17. Ionic and total abundances of NGC 6572 (PA=15°) 


Ionic abundances 

Te & Af ® 

Inner NEB 

N knot 

S knot 

Entire NEB 

He+/H+ 

1 

0.093±0.007 

0.098±0.022 

0.095±0.022 

0.096±0.007 

He + +/H + 

2 

0.0009±0.0001 

- 

- 

0.0007±0.0001 

ICF(He) 


1.00 

1.00 

1.00 

1.00 

He/H 


0.094±0.007 

0.098±0.022 

0.095±0.022 

0.097±0.007 

10^ X 0+/H + 

1 

1.52±0.11 

19.56±3.59 

8.59±1.57 

2.51±0.18 

10^ X 0^+ /H + 

2 

43.2±3.1 

34.5±5.1 

34.2±5.1 

41.9±3.1 

ICF(O) 


l.OlitO.lO 

1.00 

1.00 

l.OlitO.lO 

10^ X O/H 


4.51±0.52 

5.41±0.87 

4.28±0.66 

4.47±0.54 

10^ X N+/H + 

1 

0.67±0.06 

5.33±0.97 

3.39±0.62 

0.98±0.08 

ICF(N) 


29.6±4.2 

2.77±0.68 

4.99±1.19 

17.8±2.5 

10^ X N/H 


1.99±0.33 

1.47±0.45 

1.69±0.51 

1.75±0.28 

10 ® X S+/H + 

1 

0.13±0.02 

0.53±0.19 

0.45±0.17 

0.18±0.02 

10® X S2+ /H + 

2 

2.31±0.24 

2.92±1.05® 

2.51±0.95® 

2.34±0.24 

ICF(S) 


2.17±0.11 

1.11±0.05 

1.27±0.08 

1.85±0.08 

10® X S/H 


5.30±0.62 

3.81±1.38 

5.43±1.44 

4.65±0.52 

10® X CP+/H + 

3 

7.71±1.15 

_ 

_ 

7.67±1.07 

ICF(Cl) 


2.29±0.36 

- 

- 

1.99±0.30 

10® X Cl/H 


17.67±3.81 

- 

- 

15.24±0.31 

lo"^ X Ar®+/H + 

2 

3.77±0.39 

- 

- 

3.32±0.34 

ICF(Ar) 


1.04±0.18 

- 

- 

1.07±0.18 

lo"^ X Ar/H 


3.93±0.79 

- 

- 

3.54±0.71 

10® X Ne2+/H + 

2 

9.19±0.65 

10.03±1.84 

6.72±1.24 

9.02±0.64 

ICF(Ne) 


1.04±0.15 

1.57±0.34 

1.25±0.27 

1.07±0.15 

10® X Ne/H 


9.57±1.51 

15.71±4.48 

8.41±2.38 

9.61±1.50 


“ 1: Te[N II] and ATg [S II], 2: Tg [O III] and Afg [Ar IV] and 3: Tg [O III] and A/'g [Cl III]. 

^ It is calculated by using the equation S”t2 =g+* ^4 ^q + 2 y q^O .433 ^ from Kingsburgh & Barlow (1994). 


© 0000 RAS, MNRAS 000, [2?? 











